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ABSTRACT 

Globular clusters (GCs) are tracers of the gravitational potential of their host galaxies. Moreover, their kinematic properties may 
provide clues for understanding the formation of GC systems and their host galaxies. We use the largest set of GC velocities obtained 
so far of any elliptical galaxy to revise and extend the previous investigations (Richtler et al. 2004) of the dynamics of NGC 1399, the 
central dominant galaxy of the nearby Fornax cluster of galaxies. The GC velocities are used to study the kinematics, their relation 
with population properties, and the dark matter halo of NGC 1399. We have obtained 477 new medium-resolution spectra (of these, 
292 are spectra from 265 individual GCs, 241 of which are not in the previous data set), with the VLT F0RS2 and Gemini South 
GMOS multi-object spectrographs. We revise velocities for the old spectra and measure velocities for the new spectra, using the 
same templates to obtain an homogeneously treated data set. Our entire sample now comprises velocities for almost 700 GCs with 
projected galactocentric radii between 6 and 100 kpc. In addition, we use velocities of GCs at larger distances published by Bergond 
et al. (2007). Combining the kinematic data with wide-field photometric Washington data, we study the kinematics of the metal-poor 
and metal-rich subpopulations. We discuss in detail the velocity dispersions of subsamples and perform spherical Jeans modelling. 
The most important results are: The red GCs resemble the stellar field population of NGC 1399 in the region of overlap. The blue 
GCs behave kinematically more erratic. Both subpopulations are kinematically distinct and do not show a smooth transition. It is not 
possible to find a common dark halo which reproduces simultaneously the properties of both red and blue GCs. Some velocities of blue 
GCs are only to be explained by orbits with very large apogalactic distances, thus indicating a contamination with GCs which belong 
to the entire Fornax cluster rather than to NGC 1399. Also, stripped GCs from nearby elliptical galaxies, particularly NGC 1404, may 
contaminate the blue sample. 

We argue in favour of a scenario in which the majority of the blue cluster population has been accreted during the assembly of the 
Fornax cluster. The red cluster population shares the dynamical history of the galaxy itself. Therefore we recommend to use a dark 
halo based on the red GCs alone. 

The dark halo which fits best is marginally less massive than the halo quoted by Richtler et al. (2004). The comparison with X-ray 
analyses is satisfactory in the inner regions, but without showing evidence for a transition from a galaxy to a cluster halo, as suggested 
by X-ray work. 

Key words, galaxies: elliptical and lenticular, cD — galaxies: kinematics and dynamics — galaxies:individual: NGC 1399 



1. Introduction 

1.1. The globular cluster systems of central elliptical galaxies 

Shortly after M87 revealed its rich globular cluster system 
(GCS) (iBaumll 1 9551; [Racindl 1 9681) it became obvious that bright 
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ellipticals in general host globular clusters in mu ch larger num- 
bers than spiral galaxies dHarris & Racineiri979l) . Moreover, the 
richest GCSs are found for elhptical galaxies in the centres 
of galaxy clusters, for which M87 in the Virgo cluster and 
NGC 1399 in the Fornax cluster are the nearest example s. For 
recent reviews of the field, see iBrodie & Stradeii (l2006l) . and 
lRichtleii(l2006h . 

The galaxy cluster environment may act in different ways to 
produce these very populous GCSs. Firstly, there is the paradigm 
of giant elliptical galaxy formation by t he merging of disk galax- 
ies (e.g. iToomrd 119771 Ir enzin il l2006h . That early-type galax- 
ies can form through mergers is evident by the identification of 
merger remnants and many kinematical irregularities in ellipti- 
cal galaxies (counter-rotating cores, accreted dust and molecular 
rings). 
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In fact, starting from the bimodal colour distri bution of GCs 
in some giant ellipticals, lAshman & Zepfi d 19921) predicted the 
efficient formation of GCs in sp iral-spiral mergers, before this 
was confirmed observationally fechweizer & Seitzerlll993l) . In 
their scenario, the blue clusters are the metal-poor GCSs of 
the pre-merger components while the red (metal-rich) GCs are 
formed in the material which has been enriched in the starbursts 
accompanyin g the early merger (gaseous merger model). 

However. IForbes et al.l(ll997b pointed out that the large num- 
ber of metal-poor GCs found around giant ellipticals cannot be 
explained by the gaseous merger model. These authors proposed 
a 'multi-phase collapse model' in which the blue GCs are cre- 
ated in a pre-galactic phase along with a relatively low num- 
ber of metal-poor field stars. The majority of field stars, i.e. the 
galaxy itself and the red GCs, are then formed from the enriched 
gas in a secondary star for mation epoch. This scenario is also 
supported by the findings of lSpitler et al.l (2008) who studied an 
updated sample of 25 galaxies spanning a large range of masses, 
morphological types and environments. They confirmed that the 
spirals generally show a lower fraction of GCs normalised to 
host galaxy stellar mass than massive ellipticals, thus ruling out 
the possibility that the GCSs of massive ellipticals are formed 
through major wet mergers. Further, Spitler et al. suggest that 
the number of GC s per unit halo mass is constant - thus extend- 
ing the work by Blakeslee et al.l ( Il99 7l) who found that in the 
case of central cluster galaxies the number of GCs scales with 
the cluster mass - pointing towards an early formation of the 
GCs. 

In some scenarios (e.g. ICote et al.lfT998t iHilker et al.lfT99l 
ICote et al.l 120021: iBeaslev et al.l 12002'), the accretion of mostly 
metal-poor GCs is responsible for the richness of the GCSs of 
central giant ellipticals. 

But only during the last years it has become evident that 
accretion may be important even for GCSs of galaxies in rel- 
ativel y low density regions like the Milky Way (e.g. iHelmil 
120081) . the most convi ncing case being the GCs associated wit h 
the Sagittarius stream jlbata et al.ll2001l; iBellazzini et al.ll2003l) . 
Therefore, GC accretion should plausibly be an efficient process 
for the assembly of a GCS in the central regions of galaxy clus- 
ters. 

Given this scenario one expects huge dark matter halos 
around central giant ellipticals, perhaps even the sum of a 
galax y-size dark halo and a cluster dark halo (e.g. Ike be et al] 
Il996l) . However, dark matter studies in elliptical galaxies us- 
ing tracers other than X-rays were long hampered by the lack 
of suitable dynamical tracers. Due to the rapidly declining sur- 
face brightness profiles, measurements of stellar kinematics are 
confined to the inner regions, just marginally probing the radial 
distances at which dark matter becomes dominant. One notable 
exception i s the case of NGC 61 66, the cD galaxy in Abell 2199, 
for which iKelson et al.l ( |2002|) measured the velocity disper- 
sion profile out to a radius of 60kpc. Another one is the case 
of NGC 2974 where a n Hi disk traces the mass out to 20kpc 
dWeiimans et al.ll2008l) . 

Only with the advent of 8m-class telescopes and multi- 
object spectrographs it has become feasible to study the dy- 
namics of globular cluster systems (GCSs) of galax ies as 
distant as 20 Mpc. Earlv attempts (iHu chra & Brodie Il987t 
Grillmair et al.l Il994t ICohen & Rvzhovl ,1997: Minniti et alj 



19981 iKissler-Patig et alj Il998al) were restricted to the very 



brightest GCs, and even today there are only a handful of galax- 
ies with more than 200 GC velocities measured. Large (A^gc > 
200) s amples of GC ra dial velocities hav e been pubhshed for 
M87 (ICote et al.l l200Ti) and NGC 4472 (ICote etalJl2003h in 



Virgo and NGC 1399 in Fornax dRichtler etal.ll2004 hereafter 
Paper I). Also the GCS dynamics of the nearby (~ 4Mpc) dis- 
turbed galaxy CenA (NGC 5 128) has been studied extensively 
jPeng et al.ll2004t IWoodlev et al.ll2007h . with 340 GC velocities 
available to date. 

Investigating the kinematics and dynamics of GCSs of ellip- 
tical galaxies has a twofold objective. Firstly, kinematical infor- 
mation together with the population properties of GCs promise 
to lead to a deeper insight into the formation history of GCSs 
with their different GC subpopulations. Secondly, GCs can be 
used as dynamical tracers for the total mass of a galaxy and thus 
allow the determination of the dark matter profile, out to large 
galactocentric distances which are normally inaccessible to stud- 
ies using the integrated light. These results then can be compared 
to X-ray studies. A large number of probes is a prerequisite for 
the analysis of these dynamically hot systems. Therefore, giant 
ellipticals, known to possess extremely populous and extended 
GCSs with thousands of clusters, have been the preferred targets 
of these studies. 

Regarding the formation history of GCSs, a clear picture has 
not yet emerged. Adopting the usual bimodal description of a 
GCS by the distinction between metal-poor and metal-rich GCs, 
the kinematical properties seem to differ from galaxy to galaxy. 
For example, in M87 the blue and red GCs do not exhibit a sig- 
nificant difference in their velocity dispersion (I Cote et ani200lL 
with a sample size of 280 GCs). Together with their different 
surface density profiles. Cote et al. concluded that their orbital 
properties should be different: the metal-poor GCs have prefer- 
entially tangential orbits while the metal-rich GCs prefer more 
radial orbits. 

In NGC 4472, on the other hand, the metal-rich GCs have 
a significantly lower velocity dispersio n than their blue coun- 
terparts (found for a sample 250 GCs, ICdte et all |2003[) . and 
these authors conclude that the cluster system as a whole has 
an isotropic orbital distribution. 

NGC 1399, the object of our present study, is a galaxy, 
which, as a central cluster gal axy, is similar to M87 in many re- 
spects (see lDirsch et al.ll2003 b. Here, the metal-poor GCs show 
a distinctly higher velocity dispersion than the metal-rich GCs, 
but more or less in agreement with their different density profiles 
(with a sample size o f ~ 470 GCs). A similar behaviour has been 
found for NGC 4636 (ISchuberth et alJl2006l) . 

In most other studies, the sample sizes are still too small to 
permit stringent conclusions or the separate treatment of red and 
blue GCs, but dark matter halos have been found in almost all 
cases. 



1.2. The case of NGC 1399 

NGC 1399 has long been known to host a very populous 
globular cluster system (e.g. iDirsch et al.l 12003' and references 
therein). With the photometric study by fBassino et al. (200^ 
it became clear that the GCS of NGC 1399 extends to about 
250 kpc, which is compara ble to the core radius of the cluster 
dFerguson & Sandagelll989h . Accordingly, it has always been an 
attractive target for studying the dynamics of its GCS. One finds 
t here the larg est sample of GC velocities (469) available so far 
dRichtler et al]|2004i (Paperl). lDirsch et al.ll2004l) . It was shown 
in Paper I that blue and red GCs are kinematically different, as 
was expected from their different number density profiles: the 
red GCs exhibit a smaller velocity dispersion than the blue GCs 
in accordance with their respective density profiles. Evidence for 
strong anisotropics has not been found. The radial velocity dis- 
persion profile was found to be constant for red and blue GCs. 
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However, as we think now, this could have been a consequence 
of a velocity cut introduced to avoid outliers. A dark halo of the 
NFW type under isotropy reproduced the observations satisfac- 
torily. No rotation was detected apart from a slight signal for the 
outer blue GCs. It was shown that some of the extreme radial 
velocities in conjunction with the derived dark halo were only 
understandable if they were being caused by orbits with very 
large apogalactic distances. In this paper, we extend our investi- 
gation of the NGC 1399 GCS to larger radii (80 kpc). We simul- 
taneously revise the old velocities/spectra in order to have an ho- 
mogeneously treated sample. The case of Modifi ed Newtonian 
Dynam ics (MOND) has already been discussed in Richtl er et alj 
(120081) . where it has been shown that MOND still needs addi- 
tional dark matter of the order of the stellar mass. We do not 
come back to this issue in the present contribution. 

The GCS of NGC 1399 is very extended. One can trace the 
blue GCs out to a bout 250 kpc, the red GCs only to 140 kpc 
dBassino et al.l 12006.) . Regarding total numbers, there are only 
half as many red GCs as there are blue, suggesting that the for- 
mation of GCs in mergers is not the dominant mechanism pro- 
ducing a high specific frequency, even if in the central regions 
of a proto-cluster the merger rate is supposed to be particularly 
high. 

Following Paper I, we assume a distance modulus of 3 1 .40. 
At the distance of 19Mpc, 1" corresponds to 92 pc, and 1' cor- 
responds to 5.5 kpc. 

This paper is organised as follows: In Sect. |2] we describe 
the observations and the data reduction. The velocity data base 
is presented in Sect. [3] In Sect. H) we present the photometric 
properties and the spatial distribution of our GC velocity sample. 
The contamination by interlopers is discussed in Sect. |5] The 
properties of the line-of-sight velocity distribution are studied 
in Sect. |6] In Sect. [T] we test our GC sample for rotation. The 
line-of-sight velocity dispersion and the higher order moments 
of the velocity distributions are calculated in Sect. |8] The Jeans 
modelling and the derived mass models are described in Sects. |9] 
and[TO] The results are discussed and summarised in SectionsfTTI 
and [12] 



2. Observations and data reduction 

To study the dynamics of the GCS of NGC 1399 out to large 
radii, we use data obtained with the multi-object-spectrographs 
FORS2/MXU and GMOS-S, at the VLT and Gemini-South 
telescopes, respectively. The spectral resolution and wavelength 
coverage is similar for both datasets. Details on the mask prepa- 
ration and data reduction are given below. 

The location of the GCs on the plane of the sky is shown in 
Fig.[T] The GCs from Paper I (shown as dots) occupy the inner 
region, and open symbols represent the GCs added in the present 
study, raising the number of GC velocities to almost 700 and 
extending the radial range by almost a factor of two. 

2.1. Photometric data 

Wide-field imaging in the metallicity-sensitive Washington sys- 
tem obtained for several fields in the Fornax cluster using the 
CTIO MOSAIC camercQ forms th e basis for the pho tometry 
used in this work. The central field (iDirsch et al.ll2003]) encom- 
passes the area covered by our spectroscopic study. As in Paper I, 



' The camera is mounted in the prime focus of the 4-m Blanco tele- 
scope and the field-of-view is 36' x 36'. 
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Fig. 1. Spatial distribution of spectroscopically confirmed GCs 
with respect to NGC 1399 (0,0). Dots represent GCs from 
Paper I. Circles and squares are the new GCs, measured us- 
ing VLT/FORS 2 and Gemini/GMOS, respectively. North is to 
the top and East is to the left. The positions of NGC 1399 and 
NGC 1404 are marked by crosses. 



these data were used for target selection, and the photomet- 
ric properties of our velocity-confirmed GCs are presented in 
Sect. Ell] 

Further, the analysis of the outer fields presented by 
Bassino et al. (2006) provides the number density profiles of the 
GC subpopulations out to large radii (cf. Sect. 19.21 ) which are 
required for the dynamical modelling. 

2.2. VLT-F0RS2/MXU spectroscopy 

The spectroscopic observations of 11 masks in seven fields 
(see Table[T] for details) were carried out in visitor mode 
during three nights (2002 December 1-3) at the European 
Southern Observatory Very Large Telescope (VLT) facility on 
CeiTO Paranal, Chile. We used the FORS 2 (Focal Reducer/low 
dispersion Spectrograph) instrument equipped with the Mask 
Exchange Unit (MXU). 

Our spectroscopic target s were selected from the wide-field 
photometry by Dirsch et al. (l2003l hereafter Dh-03). The masks 
were designed using the FIMS software. To maximise the num- 
ber of objects per mask, we observed objects and sky positions 
through separate slits of 2" length. Th is strategy has prev iously 
been used for the study of NGC 1399 (iDirsch et al.l l2004 D-h04 
hereafter) and NGC 4636 (ISchuberth et al.ll2006l) . 

All observations were performed using the grism 600 B 
which provides a spectral resolution of ~ 2.5 A^ 150 km s"' (as 
measured from the line-widths of the wavelength calibration ex- 
posures). 

The reduction of the FORS 2/MXU spect r a was p erformed in 
the same way as described in IDirsch et al.l (|2004|) . so we just 
give a brief description here. After the basic reduction steps 
(bias subtraction, flat fielding, and trimming), the science and 
calibration frames were processed using the apextract pack- 
age in IRAF. The wavelength calibration was performed using 
identify. Typically 18 lines of the Hg-Cd-He arc lamp were 
used to fit the dispersion relation, and the residuals were of the 
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Table 1. Summary of VLT FORS2/MXU observations (programme ID 70.B-0174 ). The six-digit number given in the second 
column is the mask identifier as given in the image header (keyword: HIERARCH ESO INS MASK ID). The seeing values (Col. 6) 
were determined from the acquisition images taken just before the MXU exposures. The numbers given in Cols. 7 and 8 refer to the 
total number of sUts and the number of GCs found on a given mask, respectively. 



Mask 


ID 


Centre Position 
(12000) 


Exp. Time 

(sec) 


Seeing 


# Slits 


#GCs 


Night 


Mask 1 


965 705 


03:37:15.9 


-35:22:18.2 


2600 


0'.'90 


130 


9 


2002-12-01 


Mask 2 


965 813 


03:37:22.1 


-35:22:18.5 


2600 


'.'90 


132 


16 




Masks 


965 930 


03:37:34.7 


-35:31:40.9 


2600 


0'.'65 


130 


24 




Mask 4 


970 048 


03:37:40.9 


-35:31:41.9 


2600 


0'.'75 


132 


27 




Masks 


970 253 


03:38:02.8 


-35:38:15.7 


2600 


0'.'85 


131 


21 




Mask? 


971 436 


03:38:48.4 


-35:32:53.6 


2600 


O'.'SO 


128 


36 


2002-12-02 


Masks 


973 401 


03:39:21.8 


-35:21:27.2 


2600 


0'.'65 


127 


7 




Mask 9 


970 843 


03:38:33.4 


-35:40:34.1 


2600 


0'.'85 


122 


9 




Mask 10 


971 105 


03:38:39.6 


-35:40:35.3 


2600 


O'.'SO 


124 


18 




Mask 12 


970 540 


03:38:09.1 


-35:38:14.6 


2600 


1 '.'35 


129 


12 




Mask 13 


972 414 


03:38:40.5 


-35:17:24.9 


2600 


1 '.'45 


127 


10 


2002-12-03 



order 0.04 A. To perform the sky subtraction for a given GC- 
spectrum, the spectra of two or three nearby sky-slits were aver- 
aged and subsequently subtracted using the skytweak task. 

2.3. Gemini GMOS spectroscopy 

We used the Gemini Multi-Object Spectrograph (GMOS) on 
Gemini-South, and the observations were carried out in queue 
mode in November 2003 and December 2004. A total of ten 
spectroscopic masks in five fields were observed. Table|2] sum- 
marises the observations. The mask layout was defined using the 
GMOS Mask Design software. Again, we selected the GC can- 
didates from the D-i-03 photometry. We chose slits of 1" width 
and 5" length. We used the B600-i-_G5323 grating, centred on 
5500 A, giving a resolution of 0.9 A per (binned) pixel. The spec- 
tral resolution is ~ 4.5 A. 

The GMOS field of view is 5.'5 x 5.'5, and the detector array 
consists of three 2048 x 4608 CCDs aiTanged in a row (2x2 
binning results in a pixel scale of 0'.'146pixer'). Thus, with the 
chip gaps being perpendicular to the dispersion direction, two 
gaps show up in the spectra. For each mask, the observations typ- 
ically consisted of three consecutive 1300 sec exposures, which 
were bracketed by exposures of a CuAr arc lamp and screen flat 
exposures. 

The data were reduced using version 1.6 of the 
gemini.gmos IRAF-package in conjunction with a num- 
ber of customised scripts. The two prominent 'bad columns' 
on the CCD-mosaic were corrected for using fixpix, with 
the interpolation restricted to the dispersion direction. Cosmic 
ray (CR) rejection was done by combining the science expo- 
sures using gemcombine with the CR rejection option. The 
wavelength calibration was performed using gswavelength: 
Chebyshev polynomials of the 4* order were used to fit the 
dispersion relation. The number of lines used in these fits varied 
depending on the location of the slit, but typically ~ 70 lines 
were identified, and the residuals (r.m.s.) were of the order 
0.15A. We carefully inspected all calibration spectra in order 
to exclude blended lines and lines in the proximity of the chip 
gaps. In the next step, the wavelength calibration was applied 
to the science spectra using gstransform. The sky subtraction 
was done using the source-free regions of the slit, by using the 
gsskysub task in interactive mode. The final one-dimensional 
spectra were extracted using gsextract: The apertures were 



typically 1" wide, and the tracing was done using Chebyshev 
polynomials of the 4*- 8* order. 

3. The velocity data base 

In this section, we detail how we build our velocity data 
base. The radial velocities are measured using Fourier-cross - 
correlation. Coordinates, colours and magnitudes are taken from 
the D+03 Washington photometry. We anticipate here that we 
adopt C-R = 1.55 to divide blue (metal-poor) from red (metal- 
rich) GCs (cf. Sect. 14. It , and the division between foreground 
stars and bona-fide GCs is made atu - 450km s"' (cf. Sect |3.7| i. 

To obtain a homogeneous data set, we also re-measure the 
velocities for the spectra used in Paper I (the velocities are tabu- 
lated in Dh-04). 

3.1. Radial velocity measurements 

The radial velocities are obtained using the IRAF-fxcor task, 
which implements th e Fourier cross-correlation technique by 
iTonrv & Davis! (Il979h . The templates (i.e. reference spectra) are 
the FORS 2/MXU spectrum of NGC 1396 (u = 815 + 8 kms"'), 
which was already used by D+04, and the spectrum of a bright 
GC in the NGC 4636 GCSQ The latter has a heliocentric velocity 
of980+15kms ', its colour is C-/? - 1 .62 and the /^-magnitude 
is 19.9. Since the spectral resolution for both datasets is simi- 
lar, we use these templates for the FORS 2/MXU as well as the 
GMOS data. 

The cross-correlation is performed on the wavelength in- 
terval 4200 < A < 5500 A. The upper bound excludes sky- 
subtraction residuals from the most prominent telluric emission 
line at 5577 A, and the lower bound ensures that we are well 
within the region for which the FORS 2 wavelength calibration is 
reliable. For the GMOS data, we did not interpolate over the chip 
gaps. Hence these features are easily identified and excluded 
from the spectral regions used for the cross-correlation. 

Our spectral database contains velocities for 1036 spectra 
where we could identify a clear peak in the cross-correlation 
function (CCF). This number does not include obviously red- 
shifted background galaxies. 

^ Th is spectrum is part of the dataset analysed in ISchuberth et al.l 
j2009h . 
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Table 2. Summary of Gemini South GMOS observations (program IDs OS 2003B Q03 1 and OS 2004B Q029). The seeing (Col. 5) 
as estimated from the 2D-spectra of bright stars. The number of slits and the number of GCs detected on a given mask are given in 
Cols. 7 and 8, respectively. 



Mask 


Centre Position 
(J 2000) 


Exp. Time 

(sec) 


Seeing 


# Slits 


#GCs 


Night 




m2003-01 


03:39:20.2 


-35:31:58.6 


2600 


0'.'58 


31 


14 


2003-11-19 




m2003-02 


U3:39:20.z 


-35:31:58.6 


3900 


0'.'80 


21 


8 


2003-11-20 




m2003-03 


03:38:13.4 


-35:18:03.0 


3900 


0'.'73 


27 


7 


2003-11-23 


* 


m2003-04 


03:38:13.4 


-35:18:03.0 


3900 


0'.'55 


24 


3 


2003-11-23 


t 


m2004-01 


03:38:13.4 


-35:18:03.0 


2600 


0'.'66 


27 


8 


2004-12-07 


* 


m2004-02 


03:38:13.4 


-35:18:03.0 


3900 


0'.'70 


24 


4 


2004-12-07 


t 


m2004-03 


03:37:44.9 


-35:19:06.2 


3900 


0'.'57 


30 


10 


2004-12-10 




m2004-04 


03:37:44.9 


-35:19:06.2 


3900 


0'.'73 


21 


1 


2004-12-10 




m2004-05 


03:37:46.0 


-35:24:41.1 


3900 


0'.'50 


31 


14 


2004-12-11 




m2004-06 


03:37:46.0 


-35:24:41.1 


3900 


0'.'58 


25 


5 


2004-12-11 




m2004-07 


03:38:22.8 


-35:34:59.2 


3900 


0'.'58 


30 


18 


2004-12-12 




m2004-08 


03:38:22.8 


-35:34:59.2 


3900 


0'.'70 


25 


5 


2004-12-12 




(★) The masks m2003-03 and m2004-01 are identical; (t) m2003-04 and m2004-02 are identical. 
See Fig.|5]for a comparison of the velocity measurements. 



For each spectrum, we adopt as velocity the fxcor mea- 
surement with the highest value of the quality parameter !Rtd 
(w hich is inversely prop ortional to the velocity uncertainty Ay, 
see lTonrv & Davislll979i for details). For 973 spectra, both tem- 
plates yielded a velocity measurement. For the remaining 63 
spectra, only one of the templates returned a robust result. 

3.2. Velocity uncertainties of tine GCs 

The left panel of Fig.|2]shows the velocity uncertainties (as com- 
puted by fxcor) for the GCs as a function of 7?-magnitude. As 
expected, the fainter GCs have larger velocity uncertainties. Red 
and blue GCs show the same trend, yet the offset between the 
median values shows that the blue GCs, on average, have larger 
velocity uncertainties. 

The right panel shows the uncertainties versus C-R colour. 
One indeed finds that the uncertainties increase as the GCs be- 
come bluer. While this might partly be due to template mis- 
matching, the paucity of absorption features in the spectra of 
the metal-poor GCs by itself leads to larger uncertainties. We 
compared the velocity measurements of some of the bluest ob- 
jects using the spectrum of a bright blue GC as template and 
did not find any significant difference in the derived velocities 
or uncertainties compared to the results obtained with the other 
templates. 

3.3. Quality flags 

Given that our spectroscopic targets have 7?-magnitudes in the 
range 18.80 < otr < 22.97, the accuracy of the velocity determi- 
nations varies strongly, as illustrated in the left panel of Fig. |2] 
With decreasing S /N, the risk of confusing a prominent random 
peak of the CCF with the 'true' peak of the function increases. 
With the goal of weeding out spurious and probably inaccu- 
rate velocity determinations, we therefore assign quality flags 
(Class A or B) to the spectra: If at least one of the criteria listed 
below is fulfilled, the velocity measurement is regarded as 'un- 
certain' and the corresponding spectrum is flagged as 'Class B': 

- Only one template yields a velocity measurement. 

- Velocities measured with the two templates deviate by more 
than 50 km s"' 
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Fig. 2. Velocity uncertainties of the GC spectra as computed 
by fxcor. In both panels, crosses and dots represent blue and 
red GCs, respectively. Left panel: fxcor-uncertainties versus 
7?-magnitude. Narrow and wide box-plots show the data for red 
and blue GCs, respectively. Right panel: fxcor-uncertainties 
vs. C—R colour overlaid with box-plots. The long-dashed line 
at = 1.55 shows the division between blue and red GCs. 

In both panels, the short-dashed lines show the cuts used for 
assigning the quality flags (cf. Sect. l3.3l ). The boxes show the in- 
terquartile range (IQR), with the band marking the median. The 
whiskers extend to 1 .5 times the IQR or the outermost data point, 
if closer. 



- Velocity uncertainty At; > 75 km s ^ 

- Quality parameter '??TD ^ 4.5 

- Relative height of the CCF peak HGHTc 



<0.15 



- Width of the CCF peak FWHMco,t > 650 km s" 
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Fig. 3. fxcor parameters and data classification. The pan- 
els show the histograms for all (unfilled histograms) and the 
'Class A' (grey histogram bars) GC spectra. Top left: velocity 
uncertainty, cut at 75kms"'. Top right: 'Rjx), the cut is at 4.5. 
Bottom left: FWHM of the fitted cross-correlation peak, the 
cut is at 650 km s"'. Bottom right: Relative height of the cross- 
correlation peak, the cut is at 0.15. In all panels, the vertical line 
shows the corresponding cut. 



- 7?-magnitude limit: m^ > 22.6 

Figure|3] shows the distribution of the fxcor parameters for 
all GC spectra as unfilled histograms. In each sub-panel, the ver- 
tical line indicates the respective value defining the 'Class B' 
(uncertain) measurements. The 'Class A' spectra are shown as 
grey histograms. 

Assigning these quality flags to the spectra yields 723 
Class A and 313 Class B measurements. 



3.4. The new spectra 

For the new data set, velocities were determined for 477 spectra, 
179 (139) of which were obtained with CMOS, and 298 (200) 
with FORS 2, where the numbers in brackets refer to the Class A 
measurements. The slightly higher fraction of Class A spectra 
found for the GMOS data is probably due to the different treat- 
ment of the sky which, for the GMOS data, was subtracted prior 
to the extraction. 



3.5. Re-measuring the spectra from Paper I 

Our re-analysis of the spectral data set from Paper I yielded 559 
fxcor velocities (GCs and foreground stars). Our database and 
the D+04 catalogues have 503 (GCs and foreground stars) spec- 
tra in common, and the results are compared below. 



Fig. 4. Comparison to the velocity measurements by lDrrsch et al.l 
(l2004l) . Velocity difference vs. /^-magnitude. Dots and crosses 
represent red and blue GCs, respectively. The dashed (dot- 
ted) lines are drawn at +100 (+200) km s '. The squares indi- 
cate the spectra with a large (> 200km s"') deviations but a 
'Class A' velocity determination according to the criteria given 
in Sect. 13.31 and the label refers to the object 90:2 for which a 
second spectrum exists (see text for details). 



The GC catalogue by D+04 (their Table 3) Hsts the 
velocitie^ for 502 GC spectra, and the authors quote 468 as the 
number of individual GCs. This number drops to 452 after ac- 
counting for doubles overlooked in the published list. 

We determine fxcor velocities for 455 (of the 502) spec- 
tra, which belong to 415 individual GCs. With a median R- 
magnitude of 22.5, the 47 spectra for which no unambiguous 
fxcor velocity measurement could be achieved belong to the 
fainter objects in our data set. 

The difference between our measurements and the val- 
ues presented in Table 3 of D+04 are plotted against the R- 
magnitude in the upper panel of Fig. |4] For most spectra, the 
agreement is very good, but a couple of objects show dis- 
turbingly large discrepancies. For faint objects, deviations of the 
order several hundred km s ' are possibly due to multiple peaks 
in the CCF which occur in low S /A^-spectra. 

However, we also find very large (> 200, km s"') differences 
for seven 'Class A' spectra (marked with squares in Fig.|4]l. Of 
these GCs, one is also present on a second mask: object 90:2 (the 
labelled object at t; = 1688 + 31, voq4 — 817 ±27 km s"') was 
observed with GMOS (spectrum GS04-M07:171, u = 1710 ± 
34kms"'), thus confirming our new measurement. The remain- 
ing six spectra (in order of decreasing brightness: 86:19, 75:9, 
90:2, 77:84, 78:102, 81:5, 75:24, and 86:114), and the spectrum 
81:55 (for which no photometry is available) are re-classified as 
'Class B'. 

Since, for three of the very discrepant spectra (75:9, 78:102, 
and 90:2), the line-measurements {u() by D+04 lie within just 
75 km s"' of our new values, we suspect that, in some cases, ty- 
pographical errors in the published catalogue might be the cause 
of the deviations. 



^ These authors give two velocities per spectrum, is the value 
derived with fxcor, and Vf refers to the direct line measurement 



(rvidlines). 
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Table 4 in D+04 lists the velocities for 72 spectr^ of 
foreground stars with velocities in the range -350 < < 
320 km s" ' . Our data base contains measurements for 48 of these 
spectra (44 objects) . The overall agreement is good, but objects 
77:6 and 91:82 show large > 1000 kms"' deviations and are re- 
classified as GCs (spectra of Class B). For the remaining spectra, 
the velocity differences are of the order of the uncertainties. 

Finally, our database contains velocities for 56 (55 objects) 
spectra (eleven are foreground stars and 45 (44) GCs) that do not 
appear in the lists of D+04. 

The reason for this discrepancy is unknown. 

3.6. Duplicate measurements 

In this section, we use the duplicate measurements to assess the 
quality and robustness of our velocity determinations. First, we 
compare the velocities obtained when exposing the same mask 
on two different occasions. Secondly, we have objects which are 
present on more than one spectroscopic mask (but the instrument 
is the same). Then, we compare the results for objects observed 
with both FORS 2 and GMOS. Finally, we compare common ob- 
jects to values found in the literature where a different instrument 
(FLAMES) was used. 

3.6.1. Double exposures of GMOS masks: 

For the GMOS dataset, two masks (marked by an asterisk/dagger 
in Table|2]i were exposed during both observing campaigns. As 
can be seen from the upper left panel of Fig. |5] where we plot 
the velocity differences versus the /^-magnitude, the velocities 
agree fairly well within the uncertainties. The offset of lOkms ' 
is negligible and the r.m.s. is about 55 km s ' . 

3.6.2. Objects observed on different FORS 2 masks 

The upper right panel of Fig. |5] compares the velocities mea- 
sured for 82 objects present on different FORS 2 masks. The 
agreement is good, and the differences are compatible with 
the velocity uncertainties. The two outliers, GCs with devi- 
ations of the order 300 km s"' (objects 89:92=90:94, - 
21.43 mag, and 81:8=82:22, rriR = 21.0mag) are both from the 
data set analysed in Paper I. Also in Table 3 of D-i-04, the cor- 
relation velocities of these objects differ significantly (by 153 
and 369kms"\ respectively). We therefore assign these objects 
(which nominally have 'Class A' spectra) to Class B. 

3.6.3. GMOS and FORS 2 spectra 

There are 15 objects which were measured with both FORS 2 
and GMOS. All of them are GCs, and for all but two photometry 
is available. The bottom left panel of Fig. |5] shows the velocity 
differences against the /^-magnitude. The offset of x lOkms ' 
is negligible, and the r.m.s. is ^ 70kms 

3.6.4. The measurements by Bergond et al. (2007) 

Bergond et al. (l2007l B-h07 hereafter) used the FLAMES fibre- 
spectrograph on the VLT to obtain very accurate (Ay ~ 1 km s" ' ) 
velocities for 149 bright GCs in the Fornax cluster. Of these ob- 
jects, 24 (21 of which have Washington photometry) were also 



the entry for 80:99 is a duplicate of 80: 100. The object in slit 80:99 
is a GC listed in Table 3 of D+04. 



targeted in this study, and the velocities are compared in the 
bottom right panel of Fig. |5] The offset is » 15kms and 
the r.m.s. is » 40kms"'. With the exception of two outliers 
(9:71=gc216.7, 9:43=gc. 154.7), the agreement is excellent. The 
reason for the deviation of these two objects remains unknown. 

3.6.5. Accuracy and final velocities 

The repeat measurements of two GMOS masks shows that 
our results are reproducible. The absence of systematic differ- 
ences/offsets between the different spectrographs indicates that 
the instrumental effects are small. 

For the GCs for which duplicate measurements exist, we list 
as final velocity the mean of the respective Class A measure- 
ments (using the "Rtd values as weights). In case all spectra were 
classified as Class B, the weighted mean of these velocities is 
used. 

3.7. Separating GCs from foreground stars 

To separate GCs from foreground stars we plot, in the up- 
per panel of Fig. |6] colours versus heliocentric velocity. The 
data points fall into two regions: The highest concentration 
of objects is found near the systemic velocity of NGC 1399 
(1441 ± 9kms"', Paperl). These are the GCs, and we note 
that all of them have colours well within the interval used by 
D+03 to identify GC candidates (horizontal dashed lines). The 
second group of objects, galactic foreground stars, is concen- 
trated towards zero velocity and occupies a much larger colour 
range 0. The velocity histograms shown in the lower panel of 
Fig. |6] illustrate that the total sample, including those objects for 
which no photometry is available (unfilled histogram), exhibits 
the same velocity structure as the one found for the objects with 
MOSAIC photometry (grey histogram). Most importantly, the 
domains of GCs and foreground stars are separated by a gap of 
~ 100 kms '. Guided by Fig.|6l we therefore regard all objects 
within the velocity range 450 < Uheiio < 2500 km s"' as bona fide 
NGC 1399 GCs. 

3.8. Tlie finai velocity catalogue 

We determined fxcor velocities for a total of 1036 spectra. 
These objects are foreground stars and GCs. Background galax- 
ies were discarded at an earlier stage of the data analysis, and 
there were no ambiguous cases, sinc e there is a substantia l 
velocity gap behind the Fornax cluster jDrinkwater et al.|[2001h . 

Our final velocity catalogue comprises 908 unique objects, 
830 of which have MOSAIC photometry. This database con- 
tains 693 (656 with photometry) GCs, 210(174 with photom- 
etry) foreground stars and five Fornax galaxies (NGC 1404, 
NGC 1396, FCC 208, FCC 222, and FCC 1241). Of the GCs, 
47 1 have velocities classified as Class A, the remaining 222 have 
Class B measurements. The velocities for the GCs and the fore- 
ground stars are available in electronic form only. For each spec- 
trum, these tables give the coordinates, the fxcor-velocity mea- 
surement, the iRxD-parameter, and the quality flag. The final 
velocities and the cross-identifications are also given. 



^ The data set includes foreground stars which do not fulfil the colour 
criteria for GC candidates because we extracted all spectra of a given 
mask, including those from the 'positioning slits' used for the MXU 
mask alignment. 



8 



Y. Schuberth et al.: The globular cluster system of NGC 1399 



200 



100 



I ' ' ' I 
h CMOS 



> 
I 

> -100 



-200 



X Blue GCs 
• Red GCs 

O Stars 



4- 



= 




18 19 



20 21 22 

niR [mag] 



23 



18 19 



20 21 22 

mR [mag] 



23 



200 



100 



I ' ' ' I ' ' ' I 
h FORS2-GMOS 



o 
s 
o 
> 



-100 



-200 



T" 



T" 




200 



100 



s 



> -100 



-200 



.1 ' ' ' I ' ' ' I 
Bergond et al. (2007) 



T" 



T" 




18 



19 



20 21 
mR [mag] 



22 



23 



18 



19 



20 21 
mR [mag] 



22 



23 



Fig. 5. Duplicate measurements: Velocity differences vs. apparent magnitude. Crosses and dots are blue and red GCs, respectively. 
Foreground stars are shown as diamonds. The error-bars are the uncertainties of the two velocity measurements added in quadrature. 
Small symbols indicate objects where at least one of the spectra is classified as 'Class B'. In all panels, the dashed lines are drawn 
at +100 km s"'. Upper left: The GMOS observations labelled in Table|2](*, t) were carried out using the same masks. The scatter 
(r.m.s. a: 55 km s) is of the order of the velocity uncertainties, and the offset, which is ~10kms ' (dotted line), is very small, f/p/^er 
right: FORS 2 data comparison: Difference between the velocities measured for common objects on different masks. The inset 
shows the histogram of the deviations, which follow a normal distribution with cr ^ 70kms"'. The two outliers with deviations 
of X 300 km s"' are re-classified as Class B (see text for details). Bottom left: The same for objects measured using both GMOS 
and VLT. The off set is small (~ 10km s '), and the r.m.s. is a; 70km s '. Bottom right: Comparison to the FLAMES velocity 
measurements by iBergond et al.l (l2007l) . The error-bars show the uncertainties of our and their velocity measurements added in 
quadrature. The scatter is small (r.m.s. = 41 km s '), and the offset, which is -15 km s"' (dotted line), is negligible. 



4. Properties of the globular cluster sample 

4. 1 . Colour and luminosity distributions 

The luminosity distribution of our GC sample is plotted in 
the upper left panel of Fig.|2l The turn-over magnitude (TOM) 
of niR = 23.3 (D-h03) is shown for reference. It illustrates 
that our spectroscopic study only probes the bright part of 
the globular cluster luminosity function (GCLF). The distribu- 
tions fOTred_aiid_blu£_GCs^ as kernel density estimates 
(e.g. IVenables & Riplevll2002h . are very similar. The median R- 
magnitude of the GCs is = 21.75. The brightest (faintest) 
cluster has a magnitude of m^ - 18.8 (22.97). 

Can our spectroscopic dataset, which contains 656 GCs with 
known colours, be regarded as a photometrically representative 
subsample of the NGC 1399 GCS? As mentioned in Sect. |ZT] 
the colours and magnitudes of our GC sample are taken from the 
D+03 photometric study. These authors found a bimodal colour 
distribution for GCs in the magnitude range 21 < m^ < 23. The 



brightest GCs, however, were discovered to have a unimodal dis- 
tribution, peaking at an intermediate colour of C-R ^ 1.55. 

Figure|2] (upper right panel) shows the colour distribution of 
our sample. The main features described by D+03 are also found 
for the spectroscopic data set: The distribution is clearly bi- 
modal, and the blue GCs show a peak near C-R ^1.3. Following 
D+03 (and Paper I), we adopt C-R = 1.55 as the colour dividing 
blue from red GCs. 

Further, the brightest clusters (m^ < 21.1, grey histogram) 
do not seem to follow a bimodal distribution. 

Given that the brightest GCs show no signs of colour bi- 
modality, the division of these GCs into 'blue' and 'red' is be 
somewhat arbitrary/artificial. As will be shown in Sect. 15.41 the 
brightest GCs form indeed a kinematically distinct subgroup. 
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Fig. 7. NGC 1399 spectroscopic GC sample: Photometric properties and spatial distribution. Upper left: GC luminosity distribution 
(bin width - 0.2 mag). The dotted line at mu - 23.3 indicates the turn-over magnitude of the GCS. The thin solid, thick solid 
and dashed lines show the kernel density estimates for all, the blue, and the red GCs, respectively. Upper right panel: Colour 
distribution. The dashed histogram shows the distribution of all 656 velocity-confirmed GCs for which MOSAIC photometry is 
available, and the solid (unfilled) histogram are the GCs fainter than niR - 21.1, and the grey histogram shows the distribution of 
the 144 brightest (niR < 21.1) GCs. The dashed, solid and dot-dashed curves shows the respective kernel density estimate for the 
same data; a bandwidth of 0.075 mag was used (same as histogram bins). The dotted line at C-R - 1 .55 indicates the limit dividing 
blue from red GCs. 

Lower left: Radial distribution. The histogram {upper sub-panel) shows all 693 GCs with velocity measurements. The thin solid, 
thick solid and dashed lines show the kernel density estimates for all, the blue and the red GCs, respectively. The radial completeness, 
i.e. the number of GCs with velocity measurements with respect to the number of GC candidates (0.9 < C-R < 2.2) from the 
D+03 photometry is shown in the lower sub-panel. The black dots (grey squares) show the values for a faint-end magnitude limit 
of 23.0 (22.75). Lower right: Azimuthal distribution of the GCs. The position angle (PA) is measured North over East, and the 
dashed vertical lines indicate the photometric major axis of NGC 1399 (110°). The upper sub-panel shows the radial distance from 
NGC 1399 vs. PA. Crosses and dots represent blue and red GCs, respectively. The histogram {lower sub-panel) has a bin width of 
10°. The thin solid, thick solid and dashed lines show the kernel density estimates for all, the blue and the red GCs, respectively. 



4.2. Spatial distribution 

The lower left panel of Fig. [T] shows the radial distribution of 
the GCs with velocity measurements. Within the central 5'there 
are more red than blue GCs, which is a consequence of the 
steeper number-density profile of the former (see Sect. 19.21 ). At 
large radii, there are slightly more blue than red GCs. The me- 
dian (projected) distance from NGC 1399 is 5.'9, 6.'5, and 5.'4 
for all, the blue and the red GCs, respectively. The galactocen- 
tric distances of the GCs lie in the range I'.l < R < 18.'3, i.e. 
6 < ^ < lOOkpc. For comparison, the data set analysed in 
Paper! covered the range I'.l < R < 9.'2. The lower sub-panel 
plots an estimate for the radial completeness of our spectro- 
scopic sample: Dots and squares show the number of GCs with 



velocity measurements divided by the number of GC candidates 
from the D+03 photometric catalogue for two different faint-end 
magnitude limits. Considering the brighter half of our GC sam- 
ple (grey squares), the completeness lies above 50 per cent for 
radii between 2' and 8', and at about 14', it drops below 25 per 
cent. 

The bottom right panel in Fig. [T] shows the azimuthal distri- 
bution of the GCs (the position angle (PA) is measured North 
through East). The upper sub-panel plots radial distance ver- 
sus PA, and the azimuthal completeness decreases drastically 
beyond ~ 8'. The lower sub-panel shows the histogram of the 
azimuthal distribution. The paucity of GCs around -90°, 0°, 90° 
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Fig. 6. Separating GCs from foreground stars: Upper panel: C-R 
colour vs. heliocentric velocity for objects with velocities in the 
range -500 <v < 3000 km s'. The dashed lines at C-/? = 0.8 
and 2.3 show the colour interval D+03 used to define GC can- 
didates. Lower panel: The grey histogram shows the velocities 
of the objects within the above colour interval. The unfilled 
histogram also includes those objects for which no photom- 
etry is available. In both plots, the short-dashed vertical line 
at 1441 km s"' indicates the systemic velocity of NGC 1399. 
The Galactic foreground stars have a mean velocity of about 
50kms^' and exhibit a larger range in colours. The dotted 
line shows the GSR velocity vector in the direction of Fornax 
(124kms"'). The dashed line at 450 kms"' indicates the limit 
adopted to separate GCs from foreground stars. 



and 180° results from the choice of mask positions on the plane 
of the sky (cf. Fig.[T]l. 

5. Sample definition and interloper Removal 

5.1. GC subpopulations 

The principal division is the distinction between red and blue 
GCs, which, as shown in Paper I and D-i-03, respectively, behave 
differently with regard to kinematics and spatial distribution. 
However, one also has to consider that the brightest GCs seem to 
have a unimodal colour distribution (cf. Sect l4.1l D+03). To see 
how the kinematic properties of the GC populations change with 
the luminosity we plot in Fig.[8]the velocity dispersion as a func- 
tion of 7?-magnitude. We divide the spectroscopic GC sample 
into blue and red GCs and calculate the line-of-sight velocity 
dispersion as a function of 7?-magnitude using a Gaussian win- 
dow function (note that GCs brighter than = 20.0 are omitted 
because of their sparse spacing along the .jc-axis). For all three 
kernels ((Tmag - 0.15,0.2,0.3), the results are similar: For GCs 
with 20 < niK < 21, the dispersions of 'blue' and 'red' GCs 
are indistinguishable. For fainter GCs, down to about m« ^ 22, 
the blue and red GCs become well separated, and the respective 
dispersions do not appear to depend on the magnitude. For GCs 
fainter than ^ 22, however, the dispersions increase towards 
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Fig. 8. Velocity dispersion as function of /^-magnitude. For all 
GCs with niR > 20, the dispersion is calculated using Gaussian 
kernels with a width of cr^ag £ {0.15, 0.2, 0.3) mag. Dashed and 
solid curves show the values for the full and the 'Class A' sam- 
ples, respectively. The thick lines are the curves for cr^.^^ - 0.3. 
The dotted line at 22.6 mag indicates the magnitude limit for 
Class A spectra. The dashed line at 2 1 . 1 mag shows the division 
we adopt between bright and faint GCs. The ticks at the bottom 
of the panel show the /^-magnitudes of the GCs. 



fainter magnitudes - probably a result of the larger velocity un- 
certainties (see Sect. l3.2l i. The quality selection (i.e. the inclusion 
of Class B measurements, shown as dashed curves in Fig. [8]l, 
does not have any impact on the detected features. 

Guided by Fig. |8] we define the following samples for the 
dynamical analysis: 

- B(lue) and mr > 21.1 and C-R < 1.55 (256 GCs) 

- R(ed) and itir > 21.1 and C-R > 1.55 (256 GCs) 

- F(aint) and > 21.1 (512 GCs) 

- BR(right) and < 21.1 (144 GCs) 

- The full sample contains A(ll) 693 GCs with radial velocity 
measurements 

The main quantity we extract from our dataset is the line-of- 
sight velocity dispersion crjos. Since dispersion measurements 
react very severely to the presence of outliers and sampling char- 
acteristics, it is important to remove possible contaminants. 

As can already be seen from Fig. [1] our sample probably 
contains a number of GCs belonging to NGC 1404 (uheiio = 
1947 km s '). Our approach to identify these objects is detailed 
in Sect. EH 

Secondly, the presence of GCs with high relative velocities 
at large galactocentric radii (see Fig. |9l left panel) potentially 
has a large impact on the derived mass profiles. The treatment of 
these GCs with extreme velocities is discussed in Sect. 15.31 

In Sect. 15.41 we label the subsamples obtained from the sam- 
ples defined above after the interloper removal. 

5.2. NGC 1404 GC interloper removal 

In the sky, NGC 1404 is the closest (giant) neighbour of 
NGC 1399. It lies 9.'8 southeast of NGC 1399, which corre- 
sponds to a projected distance of only 54kpc. Its systemic 
velocity is 1947 km s"' (NEE0), which is inside the velocity 



* NASA/IPAC Extragalactic Database 
http : / /nedwww . ipac . cal tech . edu 
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Fig. 9. NGC 1399 GC velocities. Left panel: Velocities vs. galactocentric distance for the blue and red GC sample (shown in the up- 
per and lower sub-panel, respectively). Large and small dots show GCs with 'Class A' and 'B' velocity measurements, respectively. 
Objects brighter than R - 21.1 are marked by circles. The dashed line at 54 kpc shows the projected distance of NGC 1404. The 
solid and dashed horizontal lines indicate the systemic velocities of NGC 1399 and NGC 1404, respectively. Right: Densities for the 
data points for the blue and red samples (Bin and Riii cf. 15.41 ) shown as grey-scale and contour plot. The blue subsample appears to 
exhibit substructure which is not present in the red. 



range of GCs belonging to NGC 1399. Although NGC 1404 is 
reported to have a n unusually low specific frequency of only 
Sn - 2 d Richtler et al. 1992; Forbes et al. 1998), it is probable 
that GCs found in its vicinity belong to NGC 1404 rather than to 
NGC 1399. These GCs contaminate the NGC 1399 sample and 
would - if unaccounted for - lead us to overestimate the line-of- 
sight velocity dispersion and influence the measurement of the 
higher moments of the velocity distribution. 

In Fig. [To] (left panel) we show the radial velocities of the 
GCs versus the projected distance from NGC 1404. Within 3' 
from NGC 1404, the velocity distribution is skewed towards 
higher velocities, with two thirds of the GCs having velocities 
within 220km s"' of the NGC 1404 systemic velocity. The 
mean and median velocity of all GCs in this area are 1785 
and 1936 km s"', respectively (the corresponding values are 
indicated by the dashed and solid arrows). Beyond 3' from 
NGC 1404, the velocity field is cleai'ly dominated by NGC 1399 
GCs, and for the GCs within 3'-5' of NGC 1404, we find a mean 
(median) velocity of 1510 (1475) kms"'. Guided by these find- 
ings, we exclude all 23 (13 red and 10 blue) GCs within 3' 

16.6 kpc) of NGC 1404 from our analysis of the NGC 1399 
GCS. 

To distinguish the resulting sub-samples, we use Roman nu- 
merals as labels, assigning 'i' to the unaltered samples listed in 
Sect. 15.11 and 'ii' to the data excluding the GCs in the vicinity of 
NGC 1404. 

Beyond the scope of this paper remains the question whether 
the velocity structure in the vicinity of NGC 1404 is a supeiposi- 
tion in the plane of the sky, or if it traces a genuine int eraction be- 
tween the galaxies and their GCSs (as modelled bv iBekki et alj 
l2003h . A more complete spatial coverage of the NGC 1404 re- 
gion as well as a larger number of GC velocities in this area are 
required to uncover the presence of possible tidal structures. 



5.3. Velocity diagrams and extreme velocities 

Having dealt with the identification of GCs likely to belong to 
the neighbouring galaxy, we now study the velocity diagrams 
for the blue and red GC subpopulations. The main difference 
between red and blue GCs shown in Fig.|9]is that the red GCs are 
more concentrated towards the systemic velocity of NGC 1399, 
while the blue GCs have a larger range of velocities, implying 
a higher velocity dispersion for the latter. The red GCs occupy 
a wedge-shaped region in the diagram, suggesting a declining 
velocity dispersion. 

We further find in both sub-panels GCs with large relative 
velocities which appear to deviate from the overall velocity dis- 
tribution. The derived mass profile hinges on the treatment of 
these objects. As can be seen from Table 5 in ISchuberth et al.l 
(l2006h . the inclusion of only two interlopers at large radii is 
enough to significantly alter the parameters of the inferred dark 
matter halo. 

The right panel of Fig. |9] shows the densities of the data 
points for the samples which will be used in the dynamical anal- 
ysis. One notes that the (low-level) contours for the blue GCs 
are more irregular than those of the red GCs. Also, for radii be- 
tween 20 and 40 kpc, the blue GCs seem to avoid the systemic 
velocity of NGC 1399. The velocity structure found for the blue 
GCs appears to be much more complex than that of red GCs. 

The very concept of interlopers is inherently problematic: If 
we regard interlopers as an unbound population, we are able to 
identify only the extreme velocities and leave many undetected. 
Here we are primarily interested in a practical solution. To lo- 
cate possible interlopers, we therefore chose an approach simi- 
lar to the one described by Perea et al. ( 1990), who studied the 
effect of unbound particles on the results returned by different 
mass estimators (namely the virial mass estimator and the pro- 



12 



Y. Schuberth et al.: The globular cluster system of NGC 1399 



jected mass estimator as defined in lHeisler et al.lll985h . Perea et 
al. used a statistical 'Jacknife' technique to identify interlopers 
in their N-body simulations of galaxy clusters. Since the spatial 
distributions of the GC subpopulations around NGC 1399 do not 
follow the mass distribution (which is at variance with the un- 
derlying assumption of the virial mass estimator) , we choose the 
tracer mass estimator (TME, E vans et al.ll2003h which can be 
generalised to the case where the number density of the tracer 
population is different from the overall mass density. 

The method we use works as follows: For a set of GCs, 
we first calculate the quantity m^: 

m/v = {v^ ■ R) = Cproj ■ Mtme , (1) 

which is proportional to Mtme, the tracer mass estimate, 
and where v is the line-of-sight velocity relative to the 
NGC 1399 systemic velocity and R is the projected distance 
from NGC 1399. We then select, as potential oudier, the GC 
which has the largest contribution i.e. max((;^ ■/?), remove it from 
the list and calculate 

"^'■^ wTTjZ^^'-^- ^2) 

the expression for the remaining (A^ - 1) GCs. This procedure 
is repeated and, in the right panel of Fig. [TO] we plot the mag- 
nitude of the derivative of m^ with respect to «, the number of 
eliminated GCs. 

Obviously, selectively culling the high relative velocities at 
large distances drastically lowers the estimate for the total mass, 
resulting in large values of |dmN/dn|. Once the algorithm with 
increasing n starts removing objects from the overall velocity 
field, the mass differences between steps becomes smaller. 

As can be seen from the right panel of Fig.[TOl a convergence 
is reached after removing four GCs from the red subsample Rii. 
For the blue GCs (Bii), the situation is not as clear, and we decide 
to remove six GCs. For the bright GCs, we remove three GCs. 

As opposed to the constant velocity cuts (at Hheiio - 800 and 
2080kms ') used in Paper I, this algorithm does not introduce a 
de facto upper limit on line-of-sight velocity dispersion. 

Compared to schemes which use a jackknife to search for 
3cr deviations in a local (~ 15 neighbours) velocity field (which 
may be ill-defined in the sparsely sampled outer regions), our 
approach also works when the rejected GC is not the only de- 
viant data point in its neighbourhood. 

5.4. Defining tlie subsamples 

To assess the impact of the interloper removal discussed in the 
preceding sections, we assign the following labels: The primary 
label indicates the parent sample (All, BRight, Faint, Blue, and 
Red, cf. Sect. 15. Il l, and the Roman numerals i-iv refer to the 
sequence of interloper removal applied: 

i: Full data set, no interloper removal 

ii: GCs within 3' of NGC 1404 removed (Sect.|5^ 

m: Extreme velocities removed via TME algorithm (Sect. |5.3l l 

iv: Restriction to Class A velocity measurements 

To expand the radial coverage, we define the 'extended' 
samples Av, Bv, and Rv where we include the GC velocities 
presented by B+07 (with the Washington photometry given in 
[Schuberth et al. 2008). From the list of B+07, we use all Class A 
GCs classified as intra-cluster GCs (ICGS), plus the GCs in 
the vicinity of NGC 1399 while excluding GCs within 3'of 



NGC 1404. These GCs are shown as large triangles in the up- 
per panels of Fig. [T3] where we plot the GC velocities versus the 
distance from NGC 1399. Since most GCs in the Hst of B+07 
are brighter than mR = 2 1 . 1 , no bright-end magnitude cutoff is 
applied to the 'extended' samples. 

The basic statistical properties of the samples defined above 
are listed in Tableland discussed in the following section. 



6. The line-of-sight velocity distribution 

We expect the line-of-sight velocity distribution (LOS VD) to be 
nearly Gaussian. Deviations from Gaussianity may be caused by 
orbital anisotropics, the presence of interlopers, strong variations 
of the velocity dispersion profiles with radius, or rotation. Our 
sample, being the largest of its kind so far, allows us to test the 
statistical properties of the LOSVD in detail. 

Table[3] summarises the statistical properties of the velocity 
samples defined in Sect. 15.41 Besides the first four moments of 
the distributions (i.e. mean, dispersion, skewness, and kurtosis), 
we list the p-values returned by the Anderson-Darling (AD) and 
Kolmogorov-Smirnov (KS), and the Shapiro- Wilks (SW) tests 
for normality. With the exception of the dispersion whi ch w e 
calculate using the expressions given in Pryor & Mevlan ( 19931), 
the data listed in Table[3]were obtained using the functions of the 
el071 and nortest packages in R-statistics softwar^ for the 
higher moments and the normality tests, respectively. 

The SW test (IShapiro & Wilklll965tlRovstonlll982h is one of 
the most popular tests for normality which also works for small 
samples. The Anderson-Darling (Stephens 1974) test is a vari- 
ant of the KS-test tailored to be more sensitive in the wings of 
the distributio n. The SW and AD tests a re am ong the tests rec- 
ommended by iD'Agostino & Stephens! (Il986l) : note that these 
authors caution against the use of the KS test which has a much 
smaller statistical power. A detailed discussion of the application 
of the AD te s t in th e context of galaxy group dynamics is given 
m iHou et all ( l2009l) . These authors compare the performance of 
the AD test to the more commonly used KS andx^ tests and find 
that the AD test is the most powerful of the three. They conclude 
that it is a suitable statistical tool to detect departures from nor- 
mality, which allows the identification of dynamically complex 
systems. 

The velocity histograms for the entire, the 'bright', the red 
and the blue sample are shown in Fig.[TT] In all four sub-panels, 
we plot the corresponding samples prior to any interloper re- 
moval (unfilled histograms). The solid histogram bars show the 
data after removing GCs in the vicinity of NGC 1404 and the 
GCs identified by the MxME-algorithm (in panel a, the solid his- 
togram also excludes the bright GCs). The dashed bars show 
these samples when restricted to 'Class A' velocity measure- 
ments. The striking difference between the velocity distributions 
of red and blue GCs is that the red GCs seem to be well rep- 
resented by a Gaussian while the blue GCs appear to avoid the 
systemic velocity. Also the distribution of the bright GCs seems 
to be double-peaked. 

Below we examine whether these distributions are consistent 
with being Gaussian. 
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Fig. 10. Interloper removal. Left: Velocities vs. distance from NGC 1404. Crosses and dots represent blue and red GCs, respectively. 
Large and small symbols refer to Class A and B velocity measurements, respectively. Bright GCs {mr < 21.1 mag) are marked by 
a circle. The systemic velocities of NGC 1399 and NGC 1404 are shown as solid and dashed lines, respectively. The dotted line at 
3' is the limit we adopt to separate NGC 1404 GCs from the ones belonging to NGC 1399. The dashed and solid arrow indicate the 
mean and median velocity of the GCs within 3', respectively. The grey-scale shows the density of the data points smoothed with 
a Gaussian (cTx - 0^5, o-y - 75kms"'). Right: Interloper removal using the tracer mass estimator (TME). The magnitude of the 
derivative of rriM (cf. Eq. [T]| as function of «, the number of eliminated interloper candidates. Crosses and dots show red and blue 
GCs {R > 21.1, Duo4 > 3'), respectively. The bright GCs are shown as rectangles. The arrows indicate the number of rejected GCs 
for the different samples (Bii: 6, Rii: 4, and BRii: 3 GCs). 



6.1. Tests for normality 

Adopting p < 0.05 as criterion for rejecting the Null hypothesis 
of normality, we find that none of the bright subsamples (BRi- 
BRiv) is consistent with being drawn from a normal distribution. 

The fact that the full sample (Ai, An, Av) deviates from 
Gaussianity (at the ~ 90 per level) appears to be due to the pres- 
ence of the bright GCs: the velocity distribution of objects fainter 
than 21.1 (sample Fii is the union of Bii and Rii) cannot be dis- 
tinguished from a Gaussian. 

For the red subsamples (Ri through Riv), a Gaussian seems to 
be a valid description. Only the 'extended' sample (Rv) performs 
worse (probably because it encompasses bright GCs), where the 
hypothesis of normality is rejected at the 91% level by the KS- 
test. 

While the full blue data set and the sample after removing 
the GCs near NGC 1404 (Bi and Bii) are consistent with be- 
ing Gaussian, the outlier rejection (Bin) and restriction to the 
'Class A' velocities (Biv) lead to significantly lower p-values. In 
the case of the latter, all tests rule out a normal distribution at the 
90% level, which appears reasonable given that the distribution 
has a pronounced dip (cf. Fig.[TT] panel (c), dashed histogram). 

Below we address in more detail the deviations from 
Gaussianity, as quantified by the higher moments of the LOS VD. 



6.2. Moments of the LOSVD 

6.2.1. Mean 

Assuming that the GCs are bound to NGC 1399, the first moment 
of the LOSVD is expected to coincide with the systemic velocity 
of the galaxy, which is 1441 ± 9 km s"' (Paper I). As can be seen 
from the third column in Table|3] this is, within the uncertainties, 
the case for all subsamples, with the exception of Ri, where the 
GCs in the vicinity of NGC 1404 are responsible for increasing 
the mean velocity. 

6.2.2. Dispersion 

The second moments (Table[3] Col. 5), the dispersions, are 
calc ulated using the max imum-likelihood method presented 
by iPrvor & MevlanI (1 19931) . where the individual velocities are 
weighted by their respective uncertainties. Note that these esti- 
mates do not differ from the standard deviations given in Col. 4 
by more than 5kms"'. The interloper removal, by construc- 
tion, lowers the velocity dispersion. The largest decrease can be 
found for the blue GCs, where the difference between Bi and 
Bin is 25 km s For the red sample, the corresponding value is 
17kms"'. We only consider at this stage the full radial range of 
the sample. As will be shown later, the effect of the interloper 
removal on the dispersion calculated for individual radial bins 
can be much larger. 

The main feature is that the dispersions for the blue and 
red GCs differ significantly, with values of 333 ± 16 and 239 ± 
11 kms ' for the blue and red samples (Bm and Rm), respec- 
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Table 3. Descriptive statistics of the line-of-sight velocity distribution. Below we summarise the basic statistical properties for the 
velocity samples defined in Sect. 15.41 The samples are identified in the first column while the second column gives the corresponding 
nu mber of GC veloc i ties. C ol. 3 lists the mean velocity, and the standard deviation is given Col. 4. The dispersion crio^, as returned by 
the lPrvor & Mevlanl(ll993l) estimator, is given in Col. 5. The third and fourth moments of the velocity distributions, the skewness and 
the (reduced) kurtosis k, are listed in Cols. 6 and 7, respectively, and the uncertainties were estimated using bootstrap resampling. 
Columns 8 through 10 give the p-values returned by the Anderson-Darling, Lilliefors (Kolmogorov-Smirnov), and the Shapiro- 
Wilks (W-test) test for normality, respectively. The Roman numerals refer to the sequence of subsamples defined in Sect. 15.41 where 
(i) is the full sample, prior to interloper removal, (ii): GCs within 3' of NGC 1404 removed, (m): extreme velocities removed via 
TME algorithm (Sect. l53T i; (iv): restriction to Class A velocity measurements, and (v) includes the data from B+07. 



ID Nac 


(v) 


sd 




skew 


K 


Pad 


Pks 


Psw 


(1) (2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


No colour selection, no magnitude limit 


Ai 693 


1457±12 


312 


309+9 


-0.16+0.09 


0.12+0.15 


0.10 


0.08 


0.10 


An 670 


1446±12 


306 


302+9 


-0.18+0.09 


0.21+0.16 


0.07 


0.07 


0.08 


Av 729 


1438±11 


294 


290+8 


-0.33+0.08 


0.16+0.16 


0.01 


0.06 


0.01 


Faint (m^ < 21.1) GCs, no colour selection 


Fii 493 


1449±14 


309 


304+10 


-0.07+0.10 


0.12+0.17 


0.45 


0.46 


0.69 


Fiv 297 


1432±17 


288 


285+12 


-0.23+0.11 


-0.07+0.18 


0.28 


0.37 


0.14 


Bright (m« < 21.1) GCs 


BRi 144 


1435+25 


307 


305+18 


-0.50+0.18 


0.44+0.39 


0.05 


0.03 


0.02 


BRii 140 


1428±25 


294 


292+18 


-0.59+0.20 


0.61+0.45 


0.03 


0.03 


0.01 


BRiii 137 


1424±24 


280 


279+17 


-0.67+0.21 


0.72+0.53 


0.02 


0.04 


0.01 


BRiv 127 


1412±25 


283 


282+18 


-0.63+0.22 


0.67+0.52 


0.03 


0.02 


0.01 


Blue GCs (niR > 21.1) 


Bi 256 


1452±23 


362 


358+16 


-0.06+0.12 


-0.35+0.17 


0.65 


0.50 


0.67 


Bii 246 


1435±23 


356 


351+16 


-0.04+0.12 


-0.27+0.18 


0.72 


0.55 


0.81 


Bin 240 


1415±22 


337 


333+16 


-0.21+0.11 


-0.42+0.16 


0.22 


0.14 


0.09 


Biv 139 


1426±29 


338 


334+20 


-0.27+0.14 


-0.49+0.20 


0.10 


0.06 


0.07 


Bv 336 


1420±17 


316 


312+12 


-0.3+0.11 


-0.14+0.17 


0.18 


0.39 


0.05 


Red GCs (niR > 21.1) 


Ri 256 


1472+16 


259 


256+12 


0.02+0.14 


0.02+0.21 


0.48 


0.70 


0.70 


Rn 247 


1463+16 


254 


251+12 


0.02+0.14 


0.09+0.22 


0.37 


0.60 


0.63 


Riii 243 


1463+16 


242 


239+11 


0.04+0.13 


-0.08+0.20 


0.55 


0.80 


0.66 


Riv 158 


1437+19 


238 


234+13 


-0.02+0.14 


-0.14+0.25 


0.80 


0.63 


0.76 


Rv 329 


1453+14 


254 


251+10 


-0.2+0.15 


0.33+0.36 


0.24 


0.09 


0.12 



tively. The corresponding values given in Paper I (Table 2) are 
Cbiue - 291 ± 14 and o-fed = 255 ± 13. The agreement for the 
red GCs is good, and the discrepancy found for the blue GCs 
is due to the different ways of treating extreme velocities: In 
Paper I such objects were culled from the sample by imposing 
radially constant velocity cuts (at 800 and 2080 km s"' ), while 
the method employed here (see Sect. I5.3l l does not operate with 
fixed upper/lower limits, which leads to a larger dispersion. 

6.2.3. Skewness 

Coming back to the issue of Gaussianity, we calculate the third 
moment of the LOSVD: The skewness is a measure of the sym- 
metry of a distribution (the Gaussian, being symmetric with re- 
spect to the mean, has a skewness of zero): 

1 W r -n3 

skew = - > -!■ . (3) 

The uncertainties were estimated using a bootstrap (with 999 
resamplings). We consistently find a negative skewness for the 
bright subsamples (BRi-BRiv, Fig.fTTIfc). i.e. there are more data 
points in the low-velocity tail of the distribution. 

The blue GCs also show a negative skewness, which is sig- 
nificant for the samples Bm and Biv (cf. Fig. [TTIc). A similar 
finding was already described in Paper I. 



The distribution of the red GCs is symmetric with respect to 
the systemic velocity, as is to be expected for a Gaussian distri- 
bution (cf. Fig.Elt/). 

6.2.4. Kurtosis 

In the following, k denotes the reduced kurtosis, or kurtosis ex- 
cess (to assign the value zero to a normal distribution), i.e.: 

, N r --.4 
i \ 1 Xj — X 

K^-y - 3 . (4) 

The reduced kurtosis k for the various subsamples of our data is 
given in Col. 7 of Table[3] The uncertainties t^K were again es- 
timated using a bootstrap. The brightest GCs have (marginally 
significant) positive kurtosis values, meaning that the distribu- 
tion is more 'peaked' than a Gaussian. 

The blue samples (Bi-Biv), on the other hand, have a nega- 
tive kurtosis, indicative of a more flat-topped distribution which 
might be the result of a tangential orbital bias. 

The red GCs are consistent with k - Q, which would be ex- 
pected in the case of isotropic orbits. 

6.3. Summary of the statistical tests 

We briefly summarise the main findings from the above statisti- 
cal considerations. The velocity dispersions of the red and blue 
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Fig. 11. Velocity histograms, (a): The unfilled histogram shows 
all 693 GCs (Ai). The solid histogram bars are the faint {R > 
21.1) GCs outside the 3' radius around NGC 1404 (i.e. sam- 
ple Fn), and the dashed histogram is the same for the 'Class A' 
velocity measurements after outlier rejection (sample Fiv). (b): 
The bright {R < 21.1) GCs, BRi, BRiii, and BRiv. Panels (c) 
and (d) show the same for blue and red GCs, respectively. In 
all panels, the solid and dashed vertical lines show the systemic 
velocity of NGC 1399 and NGC 1404, respectively. The dashed 
curve is a Gaussian with the width as given in Table[3] for the 
data represented as solid histogram (i.e. Fii, BRiii, Bin, and Riii, 
respectively). Note that the range of the j/-axis has been adjusted 
for easier comparison, with panel (a) having twice the range of 
the subsequent plots. 



subsamples are, as in Paper I, significantly different, with the 
blue GCs showing a larger dispersion. We note that the disper- 
sions of the bright (R > 21.1) subsamples have values in between 
those found for the red and blue GCs. 

The entirety of the GCs (Ai,Aii) is not Gaussian, which is 
most likely due to the pronounced skewness of the brightest 
GCs. 

The blue subsamples, although formally consistent with be- 
ing Gaussian, show (after the outlier rejection) significantly neg- 
ative values for both skewness and kurtosis. 

The red GCs are well represented by a Gaussian, which 
would be expected in the case of isotropic orbits. 
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Fig. 12. Rotation of the NGC 1399 GCS. Velocities versus po- 
sition angle. The upper panel shows the data for the blue sam- 
ple Bm (i.e. the data after the removal of the GCs in the vicin- 
ity of NGC 1404 and the outlier rejection) for the azimuthally 
complete radial interval R <W.Q. Large symbols are GCs in the 
'outer' (4 < /? < 8fO) region for which a significant rotation sig- 
nal is found. The solid curve shows the best-fit sine-curve for 
this sample, the dashed curves indicate the ±2 ■ a-\os boundaries, 
where crios - 353 + 26kms"'. The vertical solid lines and ar- 
rows indicate ©o and A0o, respectively. The lower panel shows 
the same for the red GCs (sample Rm). The velocity dispersion 
of the red GCs in the 'outer' interval is o-\os - 255 ± 19 km s"' . In 
both panels, the vertical dashed lines indicate the position angle 
(photometric major axis) of NGC 1399 (110°). The respective fit 
parameters are listed in Table|4] 



7. Rotation 

The amount of rotation found in a GCS and among its subpop- 
ulations might be indicative of the host galaxy merger history. 
The presence of rotation (rather than the absence thereof, since 
one only observes projected rotation), might help constrain GCS 
formation scenarios. 

In Paper I, we did not find any signature of rotation, with the 
exception of the outer (R > 6') metal-poor GCs, for which a 
marginal rotation amplitude of 68 ± 60kms"' was detected (for 
a position angle of ©o = 140° + 39°). Our revised and extended 
data set provides the opportunity to re-examine this finding. 

We apply the method used in Paper I, which consists of fit- 
ting a sine-curve to the velocities plotted versus the position an- 
gle (as shown in Fig.[T2li: 



yrot(©) = f sys + A,ot sin(© - ©()) 



(5) 



Further, to ensure that the rotation signal is not merely an arti- 
fact of the inhomogeneous spatial and azimuthal coverage, we 
calculate the rotation amplitude for 1000 Monte Carlo (MC) re- 
alizations, in which we keep the positions of our original dataset 
but permutate the velocities. From this we determine the fraction 
/ of MC runs for which a rotation amplitude larger than the ob- 
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served one was found: i.e. the smaller the value of /, the higher 
the significance of the rotation signal. 

Table |4] lists the values for various GC subsamples. Since, 
as can be seen from Fig. |7] (bottom right panel), the azimuthal 
coverage becomes very patchy for galactocentric distances be- 
yond about 8fO (- 44 kpc), we also give the results for the radial 
range 0' < R < 8fO, where the spatial coverage is more uni- 
form and the completeness (cf. Fig.|7] bottom left panel) is rela- 
tively high. This selection is then further divided into an 'inner' 
(0' <R< 4.'0) and an 'outer' (4:0 <R< 8.'0) region. 

7.1. Rotation of the brightest GCs 

The brightest GCs < 21.1), show no sign of rotation. For all 
selections (BRi-BRiv) and radial intervals considered in Table|4] 
the rotation amplitude is consistent with being zero. 

7.2. Rotation of the metal-rich GCs 

For the metal-rich GCs, only the sample Riii (i.e. the data set 
obtained after masking the GCs in the vicinity of NGC 1404 and 
applying the outlier rejection algorithm) shows a weak rotation 
signal (A = 61 ±35kms"',0 = 154 ± 33°) for the 'outer sam- 
ple' . The data and the fit are shown in the lower panel of Fig.[T2l 
The corresponding /-value of 0.27 shows that a rotation signal 
of this magnitude arises from the randomised data with a proba- 
bility of 27 per cent, meaning no significant rotation is detected 
for the red GCs. Moreover, this weak rotation signal vanishes 
completely when only the 'Class A' velocity measurements are 
considered (Riv). 

7.3. Rotation of the metal-poor GCs 

Within the central 4!0(- 22 kpc), none of the blue GC sub- 
samples shows any rotation signal. For the 'outer samples', how- 
ever, the rotation amplitudes lie in the range 110-126kms"', 
and the values obtained for ©o agree within the uncertainties. 
For the blue GCs Bii (after removal of GCs near NGC 1404), 
the amplitude is A = 110 ± 53kms"', and ©o = 130 ± 24° 
(see Fig. [12] upper panel), which agrees with the rotation val- 
ues quoted in Paper 1. After further restricting this sample to 
'Class A' velocities, the results agree very well, although the 
smaller number of data points leads to larger uncertainties. We 
caution, however, that the probability of such an amplitude re- 
sulting from the randomised data still is about 10 per cent. 

7.4. NGC 1399 and the rotation of its GCS 

As stated above, we find no rotation for the red (metal-rich) 
GC population. In this respect, the red GCs appear to reflect 
the properties of the st ellar body o f NGC 1399 for which the 
spectroscopic study by ISagha et aTl (|2000) found only a small 
(< 30kms ') rotation along the major axis (an equally small 
signal was found for the slightly ofF-centre slit position parallel 
to the minor axis). 

For the blue GCs, a significant rotation signature is found for 
the radial range between 4 and 8 arcminutes. The axis of rota- 
tion is consistent with the photometric major axis of NGC 1399 
(R\ - 1 10°). Due to the very patchy angular coverage for radii 
beyond about 8', no statement about the rotation of the GCs be- 
yond ~ 44 kpc can be made. 

In the literature the amount of rotation is quantified in terms 
of the parameter yrot/cios, (i-S- the ratio of rotational velocity to 



velocity dispersion). For NGC 1399 we find v/cr — for the 
metal-rich GCs and u/cr < 0.3 for the metal-poor GCs. 



8. Radial velocity dispersion profiles 

The line-of-sight velocity dispersion as a function of the pro- 
jected radius is the quantity we aim to reproduce with the Jeans 
models described in Sect.|9] The histograms displayed in Fig.fTTI 
and the data given in Table[3] show that the blue and red sub- 
populations have significantly different velocity distributions, as 
is already known from Paper 1. We calculate the dispersion pro- 
files separately for both subpopulations. First, we determine the 
dispersion values using the same annular bins for both subpop- 
ulations and divide our data into radial bins covering the full 
radial range (starting at 1 .0, 3.5, 5.5, 7.5, 9.5, 12.5, 15.5, and 30.0 
arcmin). The limits of the bins are shown as dotted lines in the 
upper and middle panels of Fig. [13] 

The middle panels of Fig.[T3]show the radial dispersion pro- 
files for the red and blue GCs. Circles show the values obtained 
for the GCs fainter than niR = 21.1, prior to any interloper re- 
moval (samples Ri and Bi for red and blue GCs, respectively). 
The corresponding profiles obtained after removing the GCs in 
the vicinity of NGC 1404 and the most extreme velocities (sam- 
ples Riii and Bin) are shown as filled squares. For the extended 
samples Rv and Bv, shown as diamonds, the number of GCs in a 
bin is given in parenthesis. The dispersion profiles for the fixed 
radial bins for the samples Ri through Rv and Bi-Bv are given in 
TablelO 

Next, to obtain data points with similar statistical uncertain- 
ties, we fix the number of GCs per bin (and thereby allow for a 
larger range in radial extent of the bins). In the bottom panels of 
Fig. [13] we use a moving window containing 35 GCs to plot the 
samples Riii and Biii together with the extended samples con- 
taining the Bh-07 measurements Rv and Bv. The corresponding 
values are listed in Table lA.2l 

8.1. Dispersion of the metal-rich subpopulation 

For Rm and Rv, shown in the middle left panel of Fig.[T3] we find 
that the velocity dispersion declines for radii beyond ~ 30 kpc, 
while the values for the sample Ri (prior to any interloper re- 
jection) even show an increasing dispersion. The removal of the 
GCs in the vicinity of NGC 1404 mostly affects the 4* bin where 
the dispersion drops from 260+35 to 216 + 30 km s"' . Removing 
the three potential outliers in the 5* bin by means of the TME- 
algorithm lowers the dispersion by about 90 km s ' . 

The bottom left panel of Fig. [13] compares the moving- 
bin results for the samples Rm and Rv. Both samples show 
again a smooth decline of the dispersion for larger radii. The 
marginally higher values found for the Rv data which include 
the B+07 velocities is likely due to the inclusion of brighter 
(R < 21.1 mag) clusters which have a higher velocity disper- 
sion than the fainter red GCs. The 3"* bin for Rm has a surpris- 
ingly high dispersion of 305 ± 44 km s" ' which is 50 (67) km s" ' 
higher than the value found in the previous (next) bin. This rise 
is caused by a number of GCs at high velocities (cf. Fig. O, 
leading to a mean velocity which exceeds the systemic velocity 
of NGC 1399 by ~ 120 kms '. In the upper panel of Fig. [H 
where we show the dispersion in a semi-logarithmic plot, the 
deviation of this data point becomes even more evident. In the 
modelling presented in Sect. [TO] this data point (marked by an 
asterisk in Table |A.2l l is omitted. 
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Table 4. Rotation amplitudes and angles for the different subsamples. The first column lists the sample identifiers. Columns 2 
through 4 give the number of GCs, the parameters obtained by fitting Eq.|5]to the data over the full radial range. Column 5 gives 
the fraction / of MC runs for which Amc ^ ^rot (see text for details). Columns 6-9, 10-13, and 14-17 give the GC numbers and 
parameters obtained when restricting the radial range toO < R < 8!0, < R < 4'.0, and 4 < 7? < 8'0, respectively. 
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8.2. Dispersion oftlie metal-poor subpopulation 

For the blue GCs (shown in the right panels of Fig. [T3] ) the sit- 
uation is more complicated in so far as the dispersion profile 
is not as smooth. For the full sample Bi as shown in the mid- 
dle right panel of Fig. [13] the dispersion rises quite sharply by 
over 100 km s"' for the range between ~ 10 and 35 kpc where it 
reaches almost 400 km s"' before levelling out at ~ 360, km s"' . 
The effect of the interloper removal is strongest for the fifth bin 
where removing five GCs (four of which are in the vicinity of 
NGC 1404 and one GC is rejected by the TME algorithm, sam- 
pleBm) decreases the dispersion by more than 50 km s ^ For the 
sixth bin, however, the dispersion is again somewhat higher (al- 
though the values agree within the uncertainties). For the sample 
Bi and Bm no clear trend for the behaviour of the GCs towards 
larger radii is discernible. 

Taking into account the velocities for the outer GCs by B-i-07 
shows however, that also the blue GCs have a declining velocity 
dispersion profile (sample Bv, shown as diamonds). 

The lower right panel of Fig. [13] shows the velocity disper- 
sion profiles for the samples Bin and Bv for a moving window 
containing 35 GCs. Compared to the red GCs these profiles ap- 
pear less continuous. This behaviour is likely due to inhomo- 
geneities and substructures in the velocity field as illustrated in 
the right panel of Fig. [9] 

8.3. Comparison to tlie stellar data 

In Fig. [14] we com pare our values to the results presented by 
ISaglia et al] (I2OOOI) . These authors used Gauss-Hermite polyno- 
mials to analyse absorption-line spectra of the stellar body of 
NGC 1399 out to 97"(^ 9 kpc). The radial ranges covered by 
both data sets only have a marginal overlap. 

Stellar velocity dispersion profile: In the upper panel of Fig. [T4l 
we plot the velocity dispersion profiles for the blue and red GCs 
(Bin, Riii) obtained for moving bins of 35 GCs (same as the bot- 
tom panels in Fig.[T3]l. With the exception of the third data point, 
the red GCs appear to follow the trend shown by the stars in 
NGC 1399. Only the innermost blue data point is comparable to 



the stellar data, all subsequent bins show a much higher disper- 
sion. 

The fourth moment: In the lower panel of Fig. [14] we com- 
pare Gauss-Hermite /14 values given by Saglia et al. to the 
corresponding values for the GCs which were converted from 
the K values given in Tabl e[3] using the following relation 
dvan der Marel & Franx|[T993l) : 

K^S^-h4 (6) 

We find that the value for the blue GCs (sample Bm) lies signif- 
icantly below the stellar values, while the value for the red GCs 
(Riii) seems to agree with the h4 values of the stars. This figure 
also demonstrates that k is not additive, which is illustrated by 
the grey area indicating the value obtained when combining the 
blue (Bm) and red (Rm) sample. 

9. Jeans Models 

The mass profile of NGC 1399 is estimated by comparing the 
values derived from spherical, non-rotating Jeans models to the 
observed line-of-sight velocity dispersions. The assumption of 
spherical symmetry is justified by the near-spherical appearance 
of the galaxy. To account for the degeneracy between mass and 
orbital anisotropy, which can only be broken using much larger 
datasets than ours, we calculate models for different (constant) 
values of the anisotropy parameter /3. 

9.1. The Jeans equation 

The spherical, non-rotating Jeans equation reads: 

d{{(r)crj(r)) ^(,) „ , G ■ M(r) 

+ 2 {(r) cr^(r) = -{(r) r , (7) 

dr r 

with y8 = 1 - -V . 

0-2 

Here, r is the true 3D radial distance from the centre and 
£ is the spatial (i.e., three-dimensional) density of the GCs; cr,. 
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Fig. 13. Velocity diagrams and line-of-sight velocity dispersion profiles. Upper left: radial velocity vs. galactocentric distance for 
the red GC s. Dots and circles sh ow faint and bright {mr < 21.1 mag) GCs, respectively. Large and small triangles are the data 
points from' Bergond et aP (l2007h with Class A and B velocity measurements, respectively (see text for details). Squares mark GCs 
rejected on account of their proximity to NGC 1404; crosses are the GCs removed by the TME rejection algorithm. The solid and 
dashed horizontal lines are the systemic velocities of NGC 1399 and NGC 1404, respectively. The vertical dashed line shows the 
(projected) distance of NGC 1404. The dotted vertical lines show the radial bins for which the dispersion shown in the middle panels 
is calculated. Upper right: The same for the blue GCs. Middle panels: The dispersion profiles for fixed radial bins. Circles are the 
values obtained for the full samples (Ri and Bi) prior to any interloper removal. Filled squares show the values after removing 
GCs in the vicinity of NGC 1404 and the outliers identified by the T ME algorithm (samples Riii and Bin). The dispersion profiles 
for the extended samples (Rv and Bv) including the velocities from Bergond et al.l (|2007|) are shown as diamonds. The labels in 
parentheses give the number of GCs in a given bin. Bottom panels: Filled squares are the profiles for Riii and Bin obtained for a 
moving bin comprising 35 GCs. Diamonds show the same for the samples Rv and Bv. These moving-bin profiles which are used 
for the modelling in Sect.[10]are listed in Table lA.2l 



and cTg are the radial and azimuthal velocity dispersions, respec- 
tively, p is the anisotropy parameter, M{r) the enclosed mass and 
G is the constant of gravitation. 

For our analysis , we use the expressions given by 



e. 



Mamon & LokasI (l2005 i). see also Ivan der Marel & Franxl 
1993h . Given a mass distribution M(r), a three-dimensional 
number density of a tracer population ({r), and a constant 
anisotropy parameter/?, the solution to the Jeans equation (Eq.|7J 
reads; 

\2/J 



£{r)ai{r) 



({s)M{s)-^[^-^ ds. (8) 



This expression is then projected using the following integral: 



N{R) 



Jr 



fcr^ r dr 
Vr2 - R2 



T 

Jr 



picT^. dr 



r^Jr^-R^ 



(9) 



where A^(^) is the projected number density of the tracer pop- 
ulation, and (Tios is the line-of-sight velocity dispersion, to be 
compared to our observed values. In the following, we discuss 
the quantities required to determine crio^iR). 
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Fig. 14. Compar ison to the stellar kinematics presented by 
ISaglia et"aLl (12000). Upper panel: Line-of-sight velocity disper- 
sion profiles. Diamonds show the stellar kinematics, dots and 
crosses are our measurements of the red (Riii) and blue (Bin) 
GCs, respectively. Lower panel: Gauss-Hermite parameter /14 
(Eq.|6|l vs. projected radius. The values for the blue and red GCs 
are shown as long-dashed and dashed rectangle, respectively. 
The solid grey rectangle shows the value for the combined red 
and blue sam ple. The va l ues de rived from absorption-line spec- 
troscopy by Sagli a et akl (12000') are shown as diamonds, and the 
dashed line represents their best fit model (cf. their Fig. 5). 

Table 5. Globular cluster number density profiles: Fitting a cored 
power -law (Eq.fTOb to the data given in Table 2 of lBassino et al.l 
(I2OO6I) yields the parameters listed below (the profile for 'all' 
GCs was obtained by adding the number counts of the red and 
blue GCs). 
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Red 
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1.02±0.04 


1.98 
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8.39±1.05 
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2.32 
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35.54±6.13 


1.74±0.27 


0.84±0.02 


2.22 



9.2. Globular cluster number density profiles 

The GCS of NGC 1399 has been the target of two wide-field 
photometric studies. The work presented by (Bassino et al. 2006, 
hereafter B+06) extends the earlier D+03 study upon which the 
analysis presented in Paper I was based. For our analysis, we use 
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Fig. 15. GC number density profiles. The data from B-i-06 are 
shown together with the fits to Eq. [10] Circles and squares rep- 
resent the red and blue GCs, respectively. The fits to the red and 
the blue number density distributions are shown as dashed and 
solid line, respectively. The coefficients are listed in Table|5] The 
dashed and solid arrows indicate the core radii Rq as returned 
from the fits for the red and blue GCs, respectively. Diamonds 
show the surface brightness profile of NGC 1399, as given in 
Table 4 of D-i-03 (the data points have been shifted for compari- 
son with the red GC number density profile). 



a cored power-law profile (Reynolds-Hubble law) to fit the data 
fromBH-06: 



N(R) = A^o 



(10) 



Here, Rq is the core radius, and 2 ■ a is the slope of the power-law 
in the outer region. For the above expression, the Abel inversion 
has an analytical solu tion (i.e. Eq.[TO]can be deprojected exactly, 
e.g. lSaha et al.ll99^ . and the three-dimensional number density 
profile reads: 
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where S is the Beta function. Table|5]lists the values obtained by 
fitting Eq.[TO]to the number density profiles given in Table 2 of 
B-i-06. For the red GCs, the fit was performed for R < 35', and 
R < 45' for the blue/all GCs. 

FigurefTS] shows the number density profiles for the blue 
(open squares) and red GCs (circles) as given in B+06, together 
with the coiTesponding fits. One clearly sees that the red GCs 
have a steeper number density profile than the blue GCs. For 
reference, the surface brightness profile of NGC 1399 (as given 
by D+03, but scaled to match the red GC profile) is shown with 
diamonds. Note that in the region of overlap (2' < R < 20') the 
slope of the densities of the stars and the red GCs are indistin- 
guishable. 
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9.3. Luminous matter 

In order to assess the dark matter content of NGC 1399, we first 
need a model for the luminous matter Moreover, the spatial dis- 
tribution of the stars is needed to model the stellar velocity dis- 
persion profile (see Sect. ll0.4r2l l. 



9.3.1. Deprojection of the surface brightness profile 

Dh-03 give the following fit to the /?-band surface brightness 
profile of NGCl 399; 



fi(R) = 2.125 log 
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\o:o55/ 



+ 15.75 



(12) 



where R is the projected radius. Assuming an absolute solar lu- 
minosity of Mq j{ - 4.28, the surface brightness profile in units 
of Lopc"^ reads: 



I(R) = l.lOx lO"* 
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(13) 
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Fig. 16. Comparison between different deprojections of the sur- 
face brightness profile of NGC 1399: The dashed curve is the 
re-projection of the luminosity density profile as given in D-t-03 
(Eq.fTSTl. The solid curve shows the fit to the data (Eq.[T2li. which 
corresponds to the re-projection of Eq.[T4l 



At a distance of 19Mpc, 0.'055 correspond to 304 pc. Using 
Eq.[TT]we obtain the deprojected profile (in units of Lq pc"^): 



j(r) 



pc^ 



16.33 



1 + 



304 pc 



-1.35 



(14) 



We now compare this to the deprojected profile given in D-i-03. 
These authors obtained the three-dimensional luminosity distri- 
bution by numerically projecting a cored power-law and com- 
pare the result to the fitted surface brightness profile. The lumi- 
nosity density profile they obtained using this procedure reads: 
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(15) 



One notes that the central luminosity density of this profile is 
more than six times higher than what we find in Eq. [14] Further, 
the profile presented by D-i-03 is slightly steeper 

To show the difference between the two deprojections, we 
plot in Fig. [16] the surface brightness profile of NGC 1399 as 
given in Table 4 of D+03, together with the fit (Eq.[T2] solid line) 
and the re-projection of their luminosity density profile (dashed 
line). The latter overestimates the measured surface brightness 
of NGC 1399 and therefore results in a larger stellar mass. 

9.3.2. Stellar mass-to-light ratio: 

In Paper I, Richtler et al. quote an 7?-band stellar mass-to- 
light ratio of M/Lr - 5.5, wh ich they obtained b y convert- 
ing the B-band value derived bv'Sa glia et al.l ( l2000 f). Recently, 
[Cebhardt et al. (2007) used the Hubble Space Telescope to ob- 
tain kinematics for the very central part of NGC 1399 in order to 
measure the mass of the nuclear black hole. From their dynam- 
ical models, these authors find a best-fit M/Lr = 5.2 + 0.4 for 
a distance of 21.1 Mpc. At a distance of 19.0 Mpc, as assumed 
in this study, this corresponds to M/Lr = 5.77 ± 0.4, which is 
in excellent agreement with the value given in Paper! We thus 
adopt a stellar mass-to-light ratio of = M/Lr =5.5. 



9.3.3. Enclosed stellar mass 

We obtain the enclosed stellar mass by integrating Eq.[T4l 
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(16) 



which can be expressed in terms of the hypergeometric function: 



M(R) = C,--7:r^R^ 2F1 



R^\ 



(17) 



where C\ - 16.33 is the constant, B - 304 pc the 
core radius, and y - 1.35 is the exponent from Eq. [T4| 
2Fi(a,b;c;z) is the hypergeo metric function as defined in e.g. 
lAbramowitz & Stegunl (ll964l) . 

In practice, the hypergeometic function can be calculated 
to arbitrary precision (e.g. Press et aL 1992 ), but the computa- 
tion of the integral over M{R) becomes very time-consuming. 
In the programme calculating the Jeans-models, we therefore 
apply a piece-wise definition of the mass profile, using only 
standard functions: For radii smaller than 3430 pc, we use a 
6*-order polynomial, and a 3'^''-order polynomial in the range 
3430 < R < 9200. For R > 9200 pc the following expression is 
a good approximation to Eq. [T7] 
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(18) 



where T^o,* = 1598 pc, and C2 = 0.106. 



9.4. The anisotropy parameter /3 

The solutions to the Jeans equation given in Eq. [8] are defined 
for constant values of the anisotropy parameter. Since /3 cannot 
be determined from the data itself, we calculate our models for 
a set of anisotropy parameters. The higher moments derived in 
Sect. 18.31 suggest that the red cluster population is most likely 
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isotropic, possibly slightly radial, while the blue GCs might 
show a mild tangential bias. We therefore use [i - -0.5, 0, +0.5 
for modelling the GCs. 

9.5. Dark matter profiles 

Numerical cosmological simulation s predict c uspy dark mat- 
ter density profiles (e.g. .Bullock et al.l 1200 lb . However, the 
very inner shape appare ntly depends on numerical details 
(e.g. iDie mand et al.l l2005l) . while agreement is reached regard- 
ing the ou ter profile which declines as R'^. Here we use an NFW 
dNavarro et al.ll 19971) halo to represent a cuspy dark component. 
The density profile is given by: 



and the cumulative mass reads: 



mil. 



(19) 



(20) 



However, in low surface brightness galaxies, where dark matter 
contributes significantly already in the central regi ons, most ob- 
servations are not compatible with a cuspy halo (iGentile et al.l 
[2007) but rather with a cored density profile. We therefore also 
model a dark matter density profil e with a core. O ne possibility 
is to use the density profile, which iBurkerj (1 1 9951) introduced to 
represent the dark matter halo of dwarf galaxies. For this halo, 
the density profile is: 



Q{r) 



(21) 



and the cumulative mass is given by the following expression: 
M(r) = 











1 /Ml 
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10. Mass models for NGC 1399 



(22) 



In this Section, we present the results of the Jeans modelling. 
First, we treat the tracer populations separately, i.e. we fit the 
velocity dispersion profiles of the red and blue GC samples using 
NFW halos. In a second step, we search for a set of parameters to 
de scribe both populat ions. The stellar velocity dispersion profile 
by ISaglia etld] (l2000i) is then used to put further constraints on 
the mass model of NGC 1399. Finally, the results for the NFW 
halos are compared to cored Burkert halo models. 



tOA. Jeans analysis 

To obtain an estimate for the NGC 1399 mass profile, we com- 
pare the observed velocity dispersion profiles to Jeans models 
(Eq.|9]l. The stellar mass-to-light ratio is assumed to be constant 
(see Sect. 19.3b . and the (constant) anisotropy parameter takes the 
values /3 - -0.5, 0, or +0.5 (corresponding to a mild tangential, 
isotropic, and a slightly radial orbital bias). To find the best Jeans 
model, we adjust the parameters of the dark halo. 
The dark matter halos considered here (NFW and Burkert ha- 
los) are characterised by two parameters, a scale radius rdark 



and a density pdark- For a given tracer population and anisotropy 
/3, we calculate a grid of models where the density acts as 
free parameter while the radii have discrete values, i.e. rdark £ 
{1, 2, 3, . . . , 200) kpc. The results of this modelling are sum- 
marised in Table|6] which quotes the best-fit halo parameters and 
the corresponding values for the different tracer populations 
and anisotropics. 

The con fidenc e level (CL) contours are calculated using the 
definition bv lAvnil (Il976h . i.e. using the difference A^^ above the 
minimum value. With two free parameters, e.g. (rdark, iPdark) 
the 68, 90, and 99 per cent contours correspond to Ax^ = 
2.30,4.61, and 9.21, respectively. 

To facilitate the comparison of the results for the NFW ha- 
los to values in th e literature, we use the equations given in 
iBullock et alj (12001 ) to express the (rs,g.,) pairs in terms of the 
virial parameters (Mvh, Cyir). These authors define the virial ra- 
dius Rvit such that the mean density within this radius is Avir = 
337 times the mean density of the universe, and the concentra- 
tion parameter is defined as Cvir - /?vir/'"i- The confidence con- 
tours in Fig.[T7]are shown in the (Myir, Cyir) parameter space. 

One immediately notices that the contours for the GCs 
(shown in the top and middle row of Fig. [17) are quite elongated, 
with a negative correlation between the virial mass and concen- 
tration. Below we discuss the results obtained for the different 
tracer populations. 



10.2. NFW-models for the metal-rich population: 

The top row of Fig. [17] shows the NFW model parameters for 
the metal-rich (red) GCs. For a given value of /3, the halos de- 
rived for the sample Rm (i.e. metal-rich GCs fainter than = 
21.1, after removal of GCs in the vicinity of NGC 1404 and the 
outlier-rejection as described in Sect.O are somewhat less mas- 
sive and less concentrated than the halos found for the 'extended' 
sample Rv. The reason is that the latter includes bright GCs for 
which we found a larger velocity dispersion (e.g. Fig. [13] lower 
left panel). Note however that the models for Riii and Rv agree 
within their respective 90% contours. The results also agree well 
with the halo parameters presented in Paper I (shown as a circle 
in Fig.fTTTi. 

The modelled dispersion profiles are compared to the obser- 
vational data in the left and middle panel of Fig. [18] For both Rm 
and Rv, the agreement between models and data is very good, 
and the best-fit models (shown as thin lines) for the different 
values of the anisotropy parameter /? are indistinguishable. The 
model parameters listed in Table|6]illustrate the mass-anisotropy 
degeneracy, where the virial mass increases with /3. 



10.3. NFW models for the metal-poor GCs: 

The shape of the velocity dispersion profiles of the metal-poor 
GCs is not as smooth as that of the metal-rich GC population. 
The 'jumps' in the profiles may be caused by the presence of 
interlopers. Since the dispersion values of adjacent bins show 
quite strong variations, the x' values derived when modelling 
the blue GCs are substantially larger than those found for the red 
GCs. 

The parameters derived for the sample Bin and the 'ex- 
tended' sample Bv differ significantly: As can be seen from the 
middle row in Fig. [17] the solutions for the extended sample 
(solid contours) have a smaller virial mass and a larger concen- 
tration. 
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Table 6. Jeans modellin g best fit parameters: (a) using our GC dataset (i.e. samples Rm and Bm), (b) using t he extended data se t 
including the values from'B ergond et aLl(l2007l) (Rv and Bv), and (c) for the stellar velocity dispersion profile bv'Saglia e t al.l ( l2000l) . 
The tracer population is given in the first column, the anisotropy parameter j6 in the second. Columns 3 and 4 list the best-fit NFW 
halo parameters. The following columns give the;if^ andx^/v, where v is the number of degrees of freedom (i.e. k - 1 - p, where k 
is the number of data points and p - 2 the number of free parameters). Columns 7-9 give the virial parameters of th e dark halos, as 
defined by Bullock et al. (2001), i.e. assuming a density contrast of A = 337. The parameters for the lBurkerj (1 19951) halo are shown 
in Cols. 10-12. 
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(b) NGC 1399 GCs, extended data set including the Bergond et al. 


[2007) 


GC data (Rv and Bv): 
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(c) Stellar velocity dispersion profile 


(data from Saglia et al. 2000) 
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Compared to the models for the red GCs, the dark halos de- 
rived from the dispersion profile of the blue GCs are up to an 
order of magnitude more massive and about four times less con- 
centrated which, as will be detailed below, considerably compli- 
cates the task of finding a common solution. 

10.4. Finding joint solutions 

Given that the different tracer populations (i.e. metal-rich and 
metal-poor GCs, and stars) move in the same potential, it is nec- 
essary to find a joint solution. For a set of tracer populations 
(labelled a and b), we determine the combined parameters by 
minimising the sumx^ - xl + xl in the (rdark,edaik)-parameter 
space. Since the tracer populations may have different orbital 
anisotropics, this procedure is performed for different combina- 
tions of y6. 

1 0.4.1 . Combining red and blue GCs 

First, we assume that p - Q for both blue and red GCs. The dis- 
persion profiles for the combinations Bm + Riii (model aT) and 
the 'extended' samples Bv + Rv (model b7) are compared to the 
observations in the left and middle panel of Fig. [18] respectively: 
The thick solid lines (labelled 7B and 7R for blue and red GCs, 
respectively) do not agree with the data in the sense that the mod- 



els fail to reproduce the large difference one observes regarding 
the dispersions of red and blue GCs. 

To find a better agreement between model and data, we com- 
bine the most massive halo found for the red GCs (i.e. ySied = 
+0.5) with the least massive halo compatible with the blue GCs 
i^Hu, = -0.5). 

The modelled dispersions for blue and red GCs are shown as 
thick dashed and long-dashed lines (labelled 8B and 8R) in the 
left and middle panel of Fig. [18] Again, the agreement between 
data and model is poor. 

The bottom right panel of Fig.[T7]shows the contours for this 
combination: The joint solution (shown as a diamond) for the 
blue 'extended' sample (Bv,j8 = -0.5) and the red sample (Rv, 
P - +0.5) lies outside the 68% contour levels of the individual 
tracer populations. 

10.4.2. Red GCs and tine stellar velocity dispersion profile 

Th e stell ar velocity dispersion profile given by Saglia et 
al. (l2000h is used to put constraints on the mass profile. To be 
consistent, we use the same stellar mass profile as for the calcu- 
lations of the GC dispersion profiles (see Sect. 19. 3t . Since this 
profile is an approximation valid for R > 500 pc, we only con- 
sider the Saglia et al. data points outside 0' 15. The spatial density 
of the stars is given by Eqns.[T3]and[T4l 
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The detailed modelling by ISaglia et"an (|2000|) revealed that 
the stars have a radially variable anisotropy parameter piR), 
which is positive and almost reaches p - +0.5 near ~ 10" 
(see their Fig. 5). We therefore model the stars for the constant 
anisotropics /3 - and j3 - +0.5 (models cl and c2, respec- 
tively). The model with the radial bias gives a slightly better fit 
to the data. The stellar data cover the central ~ lOkpc, i.e. a re- 
gion with R < rs- Thus, the scale radius and hence the total mass 
of the dark halo cannot be constrained by these data. However, as 
can be seen from the bottom left panel in Fig.[T7](solid contours) 
the concentration parameter c„„ is quite tightly constrained. The 
contours for the red GCs (sample Riii, /3i-ed - 0) shown in the 
same panel (dashed contours) have a substantial overlap with 
the parameter space allowed by the stars. The contours for the 
joint solution (model alO) are shown in the bottom middle panel 

of Fig. [m 

The corresponding models are compared to the observed 
velocity dispersion profiles in the right panel of Fig. [18] The 
thick line labelled 'lOR' shows the joint solution for the red 
GCs and the thick long-dashed line shows the corresponding 
model for the stars. For stars and red GCs alike the agreement 
between the combined model (alO) and the data is excellent, and 
we adopt this solution (which also agrees very well with the re- 
sult from Paper I) as our preferred mass model for NGC 1399. 
The NFW halo has a scale parameter of r, = 34 kpc, a density 
Qs = 0.0088 Mopc"^, and the virial radius is /?vii = 5 10 kpc. 

We now assume that the cumulative mass of the dark halo 
is described by the exp ressio n given in Eq. |22] This cored halo, 
introduced by Burkert (1 19951) . has two free parameters as well, 
namely the central density go and the scale radius tq. The so- 
lutions of the Jeans equation are found in the same manner 
as outlined in Sect. 110.11 with the radii being fixed to values 
ro e {1,2, 3,,... 100) kpc. 

The best-fit parameters are given in the last columns of Tabled 
Again, the solutions found for the red GCs have substantially 
smaller scale-radii than the halos describing the blue GCs. For 
a given tracer population, the cored halo models provide a fit 
of similar quality, and we compare the derived mass profiles in 
Fig. [191 We find that for a given dispersion profile the Burkert ha- 
los are less massive than the corresponding NFW halos. Within 
the central ~ 80 kpc, which is about the radial extent of the data 
sets Riii and Bin, however, the corresponding mass profiles are 
almost indistinguishable. 

11. Discussion 

11.1. Population aspects 

The strongest evidence pointing towards the existence of 
substructure within the GCSs of elliptical galaxies is 
the bimodal distributio n of GC colours. The work s by 
e.g. iKundu & Whitmord (l2001albl) : iLarsen et al.l (1200 lb and 
iPen^ et all (|2006|) confirmed that colour bimodality is a ubiq- 
uitous feature in the GCS of giant ellipticals. As opposed to the 
Milky Way, where the GC metallicity can be directly measured 
from spectra of individual GC stars and the presence of kine- 
matically distin ct subp opulations has been firmly established 
(e.g. lZinnir985l:ICotell999l) . there ai-e but very few spectroscopic 
metallicity measurements for GCs s urrounding giant elliptical 
galaxies (e.g lCohen et al.lll998ll2003l for M87 and M49, respec- 
tively). Thus, the interpretation of the colour bimodality of el- 
liptical galaxies i n terms of a bimodal me tallicity distribution, as 
proposed by e.g. lAshman & Zepj (Il992h . has not gone unchal- 
lenged: The fundamental question in this respect is whether the 
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Fig. 20. Upper panel: Velocity dispersion as a function of C-R- 
colour. The data points show the velocity dispersion for the 
Class A GCs which are fainter than m^ - 21.1 (i.e. the union 
of samples Biv and Riv) for bins of 0.1 mag width. The labels 
show the number of GCs in a given bin. For bins containing less 
than 10 GCs, no dispersion was calculated. The dashed horizon- 
tal lines are the velocity dispersions of the samples Biv and Riv 
as given in Table[3] The solid curve shows the dispersion calcu- 
lated using a Gaussian kernel with cr = 0.1 mag. For comparison, 
the dot-dashed curve shows the same for the bright Class A GCs 
(sample BRiv). Lower panel: Velocities vs. colour Crosses and 
dots show blue and red GCs (from the samples Biv and Riv, re- 
spectively) used in the upper panel. Small squares indicate GCs 
in the too sparsely populated colour bins. In both panels, the 
dashed line at C-R - 1 .55 indicates the colour used in this work 
to divide blue from red GCs. The vertical dotted lines indicate 
the colour range for which the dispersions are calculated. 



description as two separate populations is true or whether we 
rather face a continuum of properties where the bimodal appear- 
ance is just a morphological feature in the colour dist ribution, 
caused by a non-linear colour-metallicity relation (e.g. lRichtlerl 
2006; Yoon et al. 2006). 

That in NGC 1399 blue and red GCs show a distinct dynam- 
ical behaviour was already shown in Paper I. To gain a more 
quantitative understanding of the kinematical differences be- 
tween blue and red GCs, we plot in Fig. [20] the velocity dis- 
persion as function of C-7?-colour Our sample shows a jump 
in the velocity dispersion at C-/? - 1.55 rather than a velocity 
dispersion which smoothly changes with colour. In the case of 
elliptical galaxies, this has not yet been shown before and it is 
clear evidence that we indeed face two different populations. 

It is clear however that a simple colour cut does not cleanly 
separate these two populations, but that there is some contamina- 
tion of the blue population by members of the red population and 
vice versa. Intrinsically, these populations are expected to differ 
more strongly in their kinematical and structural properties than 
our observationally motivated blue and red populations. 
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Motivated by this strong evidence for two dynamically dis- 
tinct subpopulations in the NGC 1399 GCS, we now address the 
question concerning the nature and origin of the blue (metal- 
poor) and red (metal-rich) GCs in NGC 1399. 

In the following, we discuss the idea that a major part of the 
GCS, namely the blue GCs, has its origin predominantly in the 
accretion of material, supposedly in the form of dwarf galaxies, 
during the assembly of NGC 139 9 and the e ntire Forna x clus- 
ter. T his concept is not at all new (ICdte et alj [l998; Hilk er et al.l 
11999 ) but our large sample enables us to distinguish the two sub- 
populations better than it was previously possible. 

The comparison with the Galactic GCS is instructive. The 
picture of the Milky Way's halo has undergone a dramatic 
change from an homogeneous stellar distribut ion to a st rongly 
sub-structured halo (see the review article by Helmill2008 ). The 
Sagittarius dwarf don ated six GCs to the Galactic system (e.g. 
iBellazzini et al.ll2"003h . The extragalactic origin of lj Centauri as 
a former n ucleus o f a dwarf galaxy is hardly doubted any more 
jHilker & Richtlerl 120001; iHilker et al.l [20041 iBedin et al.ll2004t 
IViUanova et al.ll2007l) . and many of the Galactic halo GCs with 
anomalous horizontal branche s (HBs) are suspected to have th eir 
origin outside the Milky Way (iLee et al ]|2007HCatelar3l2009l) . 

Given the many tidal streams detected (and yet to be de- 
tected), the hypothesis that a large part of the Galactic halo has 
fallen in has strong arguments in its favour If this applies to a 
relatively isolated spiral galaxy, we would expect the effects of 
accretion to be much more important in a giant elUptical in the 
centre of an assembling galaxy cluster. 

11.1.1. Red GCs 

The metal-rich GCs resemble the stellar field population of 
NGC 1399 under various aspects: The density profiles (shown in 
Fig. [TsT i and the velocity dispersions (Fig. [141 upper panel) are 
indistinguishable in the radially overlapping domain. Further, 
the field population and the red GCs show similar anisotropics 
in terms of I14 and k, respectively. This strongly points towards a 
common formation history. This is further supported by the spa- 
tial distribution of red GCs in elongated gala xies, which closely 
follo ws the galax y Ught, (e.g.. NG C 1380 iKissler-Patig et al.1 
[1997 : NGC 1052 Forbes et al.[|2001bh . 

Currently, the mechanisms leading to the assembly of the gi- 
ant ellipticals which now reside at the centres of galaxy clusters 
are not fully understood an d a matter of debate (e.g. Collins et al.l 
l2009HColavit¥ et al.ll2009l and references th erein). The very old 
stellar ages derived for giant ellipticals (e.g. iTrager et al.ll2000l) 
point towards their formation at high redshift followed by pas- 
sive evolution with no major episodes of star formation (or GC 
formation, see below). Such an early assembly is, for instance, 
an inherent feature of the 'revised monolithic collapse scenario' 
bv lPipino et aLl(l2008h . 

According to the available observational evidence, GC for- 
mation is enhanced in regions of high star formation rat e 



jLarsen & Richtlej2000l: IWeidner et all2004l: ISchweizedl2006l) . 
In the case of NGC 1399, the first epoch of massive cluster for- 
mation was supposedly during the first major merger or the col- 
lapse phase assembling the main body of NGC 1399. As can 
be seen from on- going mergers (e.g. IWhitmore & Schweizerl 
1 19951: IWang et al] 12004.) . the starburst commences well before 
the components are actually m erg ed (for theoretical mod els see 
e.g. iMihos & Hernquisll [19961 and Idi Matteo et al.ll2008l) . Thus, 
also in the merger scenario, one expects the formation of metal- 
rich GCs at a quite early stage. The processes leading to equi- 
librium, whether it is phase mixing or violent relaxation, then 



act equally on the stellar field population and the metal-rich 
GCs, leaving the same imprints on both populations. However, 
the same processes which lead to the observed similarities be- 
tween the field stars and the metal-rich GCs also erased all 
dynamical traces from the epoch when the GCs were formed 
- thus rendering the distinction between a monolithic collapse 
scenario and a series of (early) gas-rich mergers impossible. In 
their very deep images of NGC 1399 Tal et al. (2009) did not 
find any of the classical merger signatures (i.e. shells, isopho- 
tal twists, tidal tails, plumes and fans) - which further points 
towa rds an early assembly of NGC 1399 and its predominantly 
old (Kissler-Patig et al. 1 998bt iForbes et al.ll200Tal: iKundu et al.1 
2005;.Hempel et al...2007,) metal-rich GCs. 

11.1.2. Blue GCs 

The blue GCs are distinguished from their red counterparts by 
a shallower density profile and a higher velocity dispersion. It 
had been concluded in Paper I that the difference in the density 
profile would be enough to explain the difference in the velocity 
dispersion profile. Our present treatment, however, differs from 
the one in Paper I. Firstly, we do not impose such a strict velocity 
cut. Secondly, we use an updated density profile which is steeper 
for large radii. Because of these factors, finding a halo which 
simultaneously satisfies the constraints from the red and the blue 
GCs is more difficult. 

Moreover, the fluctuations of the velocity dispersion on small 
radial scales are larger than one would expect from Poissonian 
noise alone. The velocity uncertainties of the blue GCs are ad- 
mittedly larger but this feature persists also after imposing a 
quality selection. 

In Paper I, we identified a set of predominantly blue GCs 
with extreme velocities and showed that their apogalactic dis- 
tances must be about 500-700 kpc, i.e. significantly larger than 
the traceable extent of the NGC 1399 GCS, or the core of 
the Fornax cl uster wh ich both are of the order ~ 250 kpc 
jBassino et alji2006. and lFergusonll989l respectively). 

To explain the complex kinematical properties of the metal- 
poor GCs, we propose a scenario in which the majority of the 
blue GCs surrounding NGC 13 99 belong to this galax y's very 
extended GCS as suggested by ISchuberth et all (l2008h . In ad- 
dition, there is a population of 'vagrant' GCs which formerly 
belonged to (i) dwarf galaxies (which were disrupted in the en- 
counter with NGC 1399) or (ii) one of the more massive early- 
type galaxies in the Fornax cluster core: While the large rela- 
tive velocities of the galaxies near the cluster centre make merg- 
ers very unlikely (compared to group environments where the 
velocity dispersion is sm aller), the fly-by of a galaxy might lea d 
to tidal stripping of GCs dForbes et al.lll997l:lBekki et alj|2003h . 
Since, in general, the metal-poor GCs have a shallower num- 
ber density distribution than the metal-rich GCs, the former are 
more likely to get stripped. The observed kinematic properties 
of these stripped GCs would depend strongly on the orientation 
of the 'donor' galaxy's orbit with respect to the plane of the sky 
and, of course, the impact parameter and the effectiveness of the 
stripping. 

Assuming that the stripped GCs would form a stream which 
roughly follows the donor's orbit, GCs stripped in a galaxy- 
galaxy encounter taking place in the plane of the sky would be 
next to undetectable since their velocities would be very close to 
the systemic velocity. 

For a galaxy such as NGC 1404, on the other hand, which 
has a large relative velocity (and there are signs that its orbit is 
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very inclined, cf. iMachacek et al.ll2005h . one would expect the 
stripped GCs to have extreme velocities. 

In this picture, the observed velocity field of the blue GCs is 
the superposition of a dynamically old population and streams 
of GCs whose phase-space coordinates are still correlated with 
the orbits of the donor galaxies. Depending on the geometry, 
the observed line-of sight velocity dispersion derived from this 
composite distribution can easily exceed the intrinsic dispersion 
of the dynamically old GCs belonging to NGC 1399 itself. 

However, as pointed out in Sect. 15. 2[ a robust test of this 
scenario requires a more complete spatial coverage than pro- 
vided by our current data set in conjunction with numerical 
simulations. Assuming that the contamination by stripped GCs 
is indeed the cause for the jagged appearance of the line-of- 
sight velocity dispersion profile of the NGC 1399 blue GCs 
(cf. Fig. [13] right panels), one would expect the dispersion pro- 
files of blue GCs of ellipticals residing in lower density environ- 
ments to be more regular. NGC 4636 is such a galaxy: It is sit- 
uated in the very outskirts of the Virgo cluster (its nea rest giant 
neig hbour has a projected distance of ~ 175 kpc) and iTal et al.l 
(1200 9) found no signs for a (recent) merger Aside from two blue 
GCs with extreme velocities, the quality-selected velocity dis- 
persion profile of the metal-poor GCs declines and has a smooth 
appearance tSchuberth et al...2006„.2009[) . 

11.2. Dynamical aspects 

In the following, we discuss the stellar-dynamical quantities and 
compare the derived mass profiles to the results from X-ray 
studies and cosmological N-body simulations. 

11.2.1. Rotation 

So far, no clear picture regarding the rotational properties of the 
GCSs of giant ellipticals has emerged. For instance, M87 shows 
significant rotation for both the metal-rich and the metal-p oor 
GCs, albeit with different axes of rotation dCote et al.ll2001h . In 
NG C 4472, only the metal-poor GCs show a strong rotation sig- 
nal (ICote et al.ll2003h . For NGC 4636, no rotation is detected for 
the metal-poor GC s, while there ar e indications for a rotation 
of metal-rich GCs (ISchuberth et al.ll2b06i) (for an overview and 
discussion of the ro tational behaviour of the G CSs of several el- 
liptical galaxies see iRomanowskv et al.ll20 08'). 

In the case of NGC 1399, we find no rotation signal for the 
red (metal-rich) GCs. This is i n agreement with th e findings for 
the field population for which 'Sagli a et"aLl (l2000l) quote an up- 
per limit of Drot - 30km s '. For the blue GCs, however, the 
situation is more complicated: Although the (~ 95% CL) rota- 
tion signature found for the entire blue sample (Bi, for the full 
radial range) vanishes when culling the more uncertain (Class B) 
velocity measurements from the analysis (sample Biv), the rota- 
tion signal for the radial range 4' < R < 8' (22 < R < 44 kpc, 
i.e. the outer subsamples) appears to be robust with respect to 
interloper removal and quality selection. The values (after in- 
terloper removal; A = 110 + 53kms"',0o = 130 ± 53°) ai-e 
consistent with the rotation signature reported in Paper I. 

Unfortunately, for the outermost parts (R > 50 kpc) of 
the NGC 1399 GC S where, according to numerical simulations 
dBekki et"aLl l2005l) . the rotation amplitude of the metal-poor 
GCs is expected to be largest, the data suffer from a very incom- 
plete angular coverage, thus precluding any statement regarding 
the rotation of the GCS at very large galactocentric distances. 



1 1 .2.2. Comparison to cosmological simulations 

In Fig. [TSjwe compare our solutio ns to the results presented in 
the study of Bullock et al.l fcOOl b (cf. their Fig. 4). These au- 
thors analysed a sample of ~ 5000 simulated halos with virial 
masses in the range 1O"-1O'"^M0. The concentration parameter 
Cvir = Rvir/rdark dccrcascs with growing halo mass (solid line). 
The scatter in this relation, however, is quite large, as the dashed 
lines encompassing 68 per cent of the simulated halos show. 

We find that our best-fit combined halos with rjark - 30 kpc 
have somewhat higher concentrations than simulated halos of a 
similar mass. Yet one has to bear in mind that the numerical ex- 
periments were carried out using dark matter particles only. It is 
conceivable that the presence of baryons and dissipative effects 
might act to increase the concentration of real halos. Thus we 
conclude that our best-fit halos do not stand in stark contrast to 
the halos found in cosmological N-body simulations. 

11.3. Comparison witli X-ray mass profiles 

NGC 1399 and the Fornax cluster have, due to their proxim- 
ity and brightness, been extensively observed by several X- 
ray satellite missions, an d the d ata support the existence of 
a massive dark halo (e.g. I Jones et al. 1997; Ikebe et al. 199i 
and lPaohllo et al.l l2002). Recentlv. IChurazov et al.1 (I2008i) used 
Chandra data to demonstrate that, for the central 4^5 (^ 25 kpc), 
the gravitationa l pote ntials derived from stellar kinematics 
(cf. Saglia et al. l2000t) and X-rays agree very well. Regarding 
larger scales, the Fornax clust er is known t o be morphologically 
quite complex: For example, llkebe et al] (11996 ). using ASCA 
data, reported the discovery of substructure in the X-ray maps 
which they interpreted in terms of a galaxy-sized dark halo em- 
bedded in a cluster-wide halo. 

Us ing deep high-resolution ROSAT data, i Paolillo et al.l 
( I2OO2I) confirmed this finding and identified three different struc- 
tures: The innermost component is centred on NGC 1399 and 
dominates the inner 50"(- 4.6 kpc). Then follows a second, al- 
most spherical 'galactic halo' component, the centre of which 
lies 1' SW of NGC 1399, its radial extent is of the order 400"(- 
35 kpc). The third and most extended 'cluster' component is 
more elongated and its centre lies 5^6 northeast of NGC 1399. 

In terms of the total mass profile, the existence of different 
components leads to 'shoulder-like' features near the interface 
regions (i.e. where the gas-densities of two components become 
equal). This behaviour is clearly seen for the ra nge of mass- 
profiles derived in both studies, see Fig. 17 of iPaolillo et al.l 
(200Z) for a comparison. 

iMakishima et al.l(l2001h who analysed a sample of 20 galaxy 
clusters (including Fornax and Virgo) propose that such hierar- 
chical dark matter profiles, parametrised in terms of a central 
excess mass superimposed on a cored King-type profile, are 
in general linked to the presence of a cD galaxy. Cuspy pro- 
files, on the other hand, are found in the absence of a cD galaxy 
(e.g. Abell 1060 which has two central giant ellipticals). 

In the case of NGC 1399, the transition from the central 
component to the cluster-wide component takes place at about 
60 kpc (^ 11'), a region probed by our GC sample. 

The question arises, accordingly, which consequences such 
a 'nested' mass distribution would have on the GC dynamics. 
Would one be able to detect signs of such a halo-in-halo struc- 
ture in the velocity dispersion profile of the GCs? 

We chose the isothermal profile of Ikebe et al. (their Model 1 ) 
and, using the standard y6-model for the X-ray gas density pro- 
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Fig. 21. Co mparison to the 'nested halo' mass profile by 
llkebe et aP (| l996() . Here, we plot the dispersion profiles ex- 
pected for their Modell, see text for details. The thick curves 
are the dispersion profiles for the blue GCs, where the solid and 
dashed lines correspond to p - and /5 = -0.5, respectively. 
The thin lines show the models for the red GCs. The solid line is 
P — Q, and the long-dashed line shows the model for fi - +0.5. 
The data for the red GC samples Rm and Rv are shown as dots 
and circles, respectively. The corresponding blue samples Bin 
and Bv are shown as crosses and unfilled squares, respectively. 

files 

QgAr) = Pgas.o^l + , (23) 

calculate the total gravitating mass using 

-kT /dlno dlnT\ 

M(<r)= - + \r (24) 

GyumH \ d In r d In r / 

dFabricant et al.|[T980l) . Here,ySx is the power-law exponent, 

Mh is the mass of the hydrogen atom, and fi - 0.6 the mean 
molecular weight (assuming a fully ionised plasma with primor- 
dial element abundances). The total density is the sum of two 
components with the parameters given in Table 1 of Ikebe et al.: 
£)i,o = 2.3 X IQ-^cmr^, Q2fl = 8.2 x lO^^^cm^^ r,,i = 4.8 kpc 
rc.2 - 127kpc, y6x,i - 0.51, and fix,! - 0.60. Since we assume 
that the gas is isothermal {T — 1.2 x 10^ K), the second term 
in Eg. l24l vanishes. Figure|2T| shows the velocity dispersion pro- 
files expected for the blue and red GC subpopulations (assum- 
ing constant anisotropy parameters of -0.5, 0, +0.5) for the mass 
profile given in Eq.|24]and plotted in Fig.|22] Indeed, the mod- 
elled velocity dispersion profiles show a dip near ~ 40 kpc, but 
the feature is quite shallow and appears to be at variance with 

the findings for the red GCs. 

The mass profiles derived by iPaolillo et al.l (|2002|) would 
lead to more pronounced features in the predicted velocity dis- 
persion profiles. Note however that a simple isothermal profile 
based on the gas density distribution ignoring the offsets of the 
different components (shown as thick line in their Fig. 15) leads 
to a non-monotonic (i.e. unphysical) density-profile at the tran- 
sition from the galaxy to the cluster component. 

Figure|22] compares the mass profiles derived in this study 
(shown as solid and dot-dashed lines) with the X-ray mass pro- 
files presented by Paolillo et al. (2002) and Ikebe et al. (1996). 
Within the central 100 Mpc, where the vast majority of our dy- 
namical probes is found, the preferred mass profiles derived in 
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Fig. 22. Comparison to X-ray measurements. The grey area 
shows the range of the mass models presented by iPaolillo et al.l 
(12002^ Th e thic k solid line shows the mass derived by 
Ike be et al.l (Il996h (their Model 1). The thin solid line indicates 
the mass profile derived in Paper I. The dash-dotted curves are 
the two halos presented in Ri chtler et a l. (2008): The more mas- 
sive one (labelled Rl) has the parameters - 50 kpc and 
Qs - 0.0085 Mq pc~ ^. The halo derived for the 'safe' sam- 
ple in iRichtleret all (|2008) (R2) has r, = 50 kpc and q, = 
0.0065 Mopc "^. The dashed and long-dashed curves are the 
best-fit NFW mass profiles alO and ^770 derived in this work 
(cf. Table|6]l. The stellar mass is shown as a short-dashed line. 

this study agree well with the 'I kebe et al.l (1996') mass model, 
although we do not find any signs for the transition from galaxy 
to cluster halo in our kinematic data. 

12. Conclusions 

Using the largest data set of globular cluster (GC) radial 
velocities available to date, we revisit and extend the inves- 
tigation of the kinematical and dynamical properties of the 
NGC 1399 GC system with respect to Richtler et al. ( |2004 
Paper I). We measure about 700 GC radial velocities out to ap- 
proximately 14.5' ^ 80 kpc (in Paper I we reached 40 kpc). To 
this sample we add 56 GC velocities from Bergond et al. (200^ 
B+07) which go as far out as 200 kpc. 

Our main findings are the following: There is no significant ro- 
tation signal among the red (metal-rich) subpopulation. We find 
rotation around the minor axis for the blue (metal-poor) clusters 
in a radial interval of 4' < < 8' (i.e. 22 kpc < R < 44 kpc), 
however weak. The blue and red clusters form two kinemati- 
cally distinct subpopulations rather than showing a continuum 
of kinematical properties. 

The red clusters correspond under various kinematical as- 
pects to the stellar field populations. 

Their velocity dispersion declines outwards and their 
velocity distribution suggests orbital isotropy. 

The jump to the higher dispersion of the blue clusters occurs 
at a colour (C-R - 1.55), which is also suggested by the mor- 
phology of the colour distribution. The blue GCs, however, show 
a more complex behaviour. Their velocity dispersion profile is 
not as smooth as that of the red GCs. There exist very high/low 
individual velocities, which suggest very large apogalactic dis- 
tances, as already found in Paper I. These objects seem to be- 
long to an intergalactic cluster population, which may be made 
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up of GCs stemming from disrupted dwarfs and GCs stripped off 
neighbouring galaxies during close encounters. 

We performed a Jeans analysis in order to constrain the mass 
profile (stellar plus dark). We found it difficult to find a dark halo 
which simultaneously accounts for the red and the blue clusters, 
which again argues for a different dynamical history. 

The dark NFW-halo which was found to represent the kine- 
matics of the red GCs (Pcc s = 0) and t he stellar velocity dis- 
persion profile presented bv lSaglia et alj ([2000) (Pstm-s - +0.5), 
our Model aiO, has the parameters pdark - 0.0088 Mopc"^ and 
'"dark - 34kpc. The virial mass is Mvir - 8.0 x IO^^Mq. 
Including the velocities from B+07, {Model blO) yields a 
slightly more massive halo (gdark - 0.0077 Mq pc^^ and rdark - 
38 kpc, and a virial mass Mvir = 9.5 X lO^^Mo). 

These NFW halos, which were found to represent the obser- 
vations, are marginally less massive than the one from Paper I. 
The total mass profile fits agree reasonably well with the X- 
ray based mass profiles out to 80 kpc. However, our model halos 
when extrapolated stay significantly below the X-ray masses. 
Moreover, we do not see the transition from the galaxy halo to 
the cluster halo claimed to be present in the X-ray data. 

We argue that these findings are consistent with a sce- 
nario where the red GCs are formed together with the bulk 
of the field population, most probably in early multiple merg- 
ers with many progenitors involved. This is consistent with the 
NGC 1399 GCs bei ng pred ominantl y old (Kissler-Patig et al. 
1998b; Forbes et alJ 12001 at iKundu etaP 120051: iHempel et all 
2007.) . with ag es very similar to the stellar age of NGC 1399 
(1 1.5 + 2.4Gvr. lTrager et all2000h 

A large part of the blue GC population has been acquired by 
accretion processes, most plausibly through dwarf galaxies dur- 
ing the assembly of the Fornax cluster. Moreover, there should 
be a significant population of intra-cluster GCs. 
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Table A.l. Velocity dispersion profiles for fixed annular bins (cf. Fig. [13] middle panels). The bins start at R = 
1.'0,3:5,5:5,7.'5,9.'5, 12.'5, 15.'5, and 30.'0; at a distance of 19Mpc, l.'O corresponds to 5.2 kpc. The left and right table show the 
values obtained for the blue and red GCs, respectively. The Roman numerals refer to the sample definition in Sect. 15.41 Column (1) 
gives the bin number, the mean radius of the GCs in a bin is given in Col. (2). The velocity dispersion (in units of kms"') and its 
uncertainty are given in Cols. (3) and (4). Col. (5) is the number of GCs in a given bin, and Col. (6) their mean velocity. 
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Fig. 17. Mass models (NFW halo): Confidence level contours (68%, 90%, and 99%). The jc-axis plots the virial mass (M„,>, in units 
of Mq) and the (/-axis is the dimensionless concentration parameter c„,>- Top row: Red GCs: From left to right, the panels display the 
results for/3 - -0.5, 0, and +0.5. Dashed contours show the values obtained for our GC data (Rm, models al through a3), while the 
solid contours show the extended sample (Rv, models bl through b3). The best-fit parameters for Riii and Rv are marked by a square 
and a cross, respectively. Middle row: The same for the blue GCs: samples Bm (models a4-a6) and Bv (models b4-b6). Bottom 
left: Parameters of the NFW halo calculated from the stellar velocity dispersion of NGC 1399 as published by Saglia et al. (200(J, 
assuming fi - +0.5 (solid contours, model c2) compared to the red GCs Rm for y6 = (dashed, model a2). The bottom middle panel 
shows model alO, the joint solution for the isotropic fi -Q red (Rm) and fi - +0.5 stellar model shown in the bottom left panel. The 
bottom right panel illustrates the difficulty to find a common solution for red and blue GCs. Solid lines are the contours for fi = 0.5 
for red extended sample (Rv centred on the dot) and for the blue GCs assuming fi - -0.5 (sample Bv, centred on the cross). The 
diamond indicates the joint solution, i.e. the minimum x]ed ^liue' corresponding contours are shown as dashed lines. Note 
that the best-fit joint solution lies outside the respective 68% CL contours of the individual GC tracer populations. 
In all p anels, the long-dashed (dot-dashed) lines give the median (68 per cent values) for simulated halos as found bv lBullock et al.l 
(1200 ll) (cf. their Fig. 4). For comparison, the circle (Mvh = 8.9 x 10'^ Mq, Cvh - 15.9) indicates the parameters of the best-fit NFW 
halo in Paper I. Note that the wiggles and gaps in the contours are an artifact of the finite grid-size and the re-gridding onto the 
(M„,>, c„,>)-plane. 




Fig. 18. Best fit Jeans models (NFW halo). Left panel: NGC 1399 GC data and modelled velocity dispersion profiles. Crosses 
and dots show the red (Riii) and blue (Bin) velocity dispersion profiles (for a moving bin of 35 GCs, cf. Fig. [13] lower panels), 
respectively. The thin curves are the best fit models obtained for the red (models al-a3) and the blue GCs (models a4-a6) separately. 
The thick lines (with labels) show the dispersions for the combined model for red and blue GCs (models a7 and a8). Middle panel: 
The same for the 'extended' samples including the Bergond et al. (2007) measurements (samples Rv and Bv). Right panel: Red GCs 
and stars. Diamonds show the stellar velocity dispersion (iSaglia et al.ll2000l) . the red GCs (Rm) are shown as dots. The thin lines are 
the models for the individual tracer populations, i.e. models cl and c2 for the stars and models al-a3 for the red GCs. The thick 
curves are the joint solutions, i.e. models a9 and alO. In all panels, solid lines are isotropic models, and the models with fi - -0.5 
and p - +0.5 are shown as dashed and long-dashed lines, respectively. The parameters of the halos shown here are listed in Table|6] 
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Fig. 19. Left panel: Mass profiles derived for the GC dispersion profiles Riii and Bm, i.e.. The thick lines with labels referring to 
the identifiers in Table|6]show the mass assuming an NFW-type dark halo. The cored Burkert halos are shown as thin curves with 
the same line style. For a given tracer population, the Burkert halo is less massive than the best-fit NFW halo. Right panel: The 
same for the extended data sets Rv and Bv (i.e. including velocities from B+07). In both panels, the dotted and dashed vertical lines 
indicate the radial distance of the outermost velocity dispersion data point for the red and blue GCs, respectively (cf. Table |A.2l l. 
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Table A.2. Velocity dispersion profiles for the subsets of our data used for the Jeans modelling. The dispersion profiles are obtained 
using a moving bin of 35 GCs. The bins do not overlap, and the minimum number of GCs per bin was set to 20, hence the large 
number of GCs in the last bin of sample Rv. For each bin, R is the mean radial distance from NGC 1399, and /?niin and /?niax give 
the radial range covered by the GCs of the bin (at a distance of 19 Mpc, I'O corresponds to 5.52 kpc). A^gc is the n umber of GCs in 
a give n bin. The line-of-sight velocity dispersion cr and its uncertainty Act are calculated using the estimator by iPrvor & MevlanI 
(Il993h . The asterisk marks the data point which was omitted when modelling the sample Riii. 
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